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Abstract 

We consider a Hamiltonian H = H°{p) + KH^{p,q,t), {p,q) G R" x T", t G K where k G R is a 
small perturbation parameter and p, q are the action and angle variables respectively. The Hamiltonian 
generates an autonomous vector field obtained by extending the phase space making t a dependent variable 
and adding its conjugate variable r. In this paper we look at a time aperiodic perturbation H^{p,q,t) 
which tends as t ^ oo to either a time independent perturbation or a time quasiperiodic perturbation 
and we prove a KAM-type theorem. Extending the phase space results in the preservation under a small 
enough perturbation of cylinders of the extended autonomous system rather than the usual tori. To 
prove the theorem we transform the Hamiltonian H to a normal form which depends on fewer angles, 
none if possible. This transformation is done via a near identity canonical transformation. The canonical 
transformation is constructed using the Lie series formalism and by solving for a generating function. 
Because of the aperiodic time dependence, the usual Fourier series methods used to obtain the generating 
function no longer apply. Instead, we use Fourier transform methods to solve for the generating function 
and make use of an isoenergetic non-degeneracy condition which results in a shift of frequencies associated 
with each cylinder. 



*This research was supported by ONR Grant No. N00014-97- 1-0071. 
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1 Introduction 



Our motivation for this work is the apphcation to the existence of "flow barriers" in two dimensional, 
incompressible, time-dependent fluid flows. Over the past 10 years there has been much work in applying the 
approach and methods of dynamical systems theory to the study of transport in fluids from the Lagrangian 
point of view. Suppose one is interested in the motion of a passive tracer in a fluid (e.g. dye, temperature, 
or any material that can be considered as having negligible effect on the flow), then, neglecting molecular 
diffusion, the passive tracer follows fluid particle trajectories which are solutions of 



x = v{x,t). (1.1) 
where v{x, t) is the velocity field of the fluid flow, x e R", n = 2 or 3. When viewed from the point of view 



of dynamical systems theory, the phase space of (1.1) is actually the physical space in which the fluid flow 



takes place. Evidently, "structures" in the phase space of (1.1) should have some influence on the transport 
and mixing properties of the fluid. Babiano Q and Aref and El Naschie provide recent reviews of this 
approach. 

To make the connection with the large body of literature on dynamical systems theory more concrete let 
us consider a less general fluid mechanical setting. Suppose that the fluid is two-dimensional, incompressible, 
and inviscid. Then we know that the velocity field can be obtained from the derivatives of a scalar valued 
function, tp{xi, X2,t), known as the streamfunction, as follows 



Xl = -7\ {Xi,X2,t), 

0x2 

X2 = -p^{xuX2,t), ixi,X2)eR''. (1.2) 

0x1 



In the context of dynamical systems theory, (1.2) is a time-dependent Hamiltonian vector field where the 



streamfunction plays the role of the Hamiltonian function. If the flow is time-periodic then the study of (1.2) 
is typically reduced to the study of a two-dimensional area preserving Poincare map. Practically speaking, 
the reduction to a Poincare map means that rather than viewing a particle trajectory as a curve in continuous 
time, one views the trajectory only at discrete intervals of time, where the interval of time is the period of 
the velocity field. The value of making this analogy with Hamiltonian dynamical systems lies in the fact 
that a variety of techniques in this area have immediate applications to, and implications for, transport 
and mixing processes in fluid mechanics. For example, the persistence of invariant curves in the Poincare 
map (KAM curves) gives rise to barriers to transport, chaos and Smale horseshoes provide mechanisms for 
the "randomization" of fluid particle trajectories, an analytical technique, Melnikov's method, allows one 
to estimate fluxes as well as describe the parameter regimes where chaotic fluid particle motions occur, a 
relatively new technique, lobe dynamics, enables one to efficiently compute transport between qualitatively 
different flow regimes. 

Extending the constructions of Smale horseshoes, Melnikov's method, and stable and unstable manifolds 
of hyperbolic trajectories to vector fields with aperiodic time-dependence has been done. However, extensions 
of KAM theory to vector fields with arbitrary time-dependence has not been done. This paper represents a 
first step in that direction. 

One of the main stability results for nearly-integrable Hamiltonian systems is Kolmogorov's theorem 
concerning the preservation of a set of fuU-measure-nonresonant invariant tori. Such nearly-integrable 
Hamiltonian systems are generated by a real valued Hamiltonian written in action-angle variables 

H[p,q) = H\p) + KH\p,q), 

where p = (pi, ...,p„) e B C M", q = (gi, ...,g„) £ T" are, respectively, the action and angle variables, and 
K is a small perturbation parameter. KH^(jj,q) is therefore a small perturbation of the integrable Hamil- 
tonian H'^{p). For any p G B the unperturbed angular frequencies are defined by X{p) = (Ai, An)(p) = 
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{dH^ /dpi{p), ...,dH^ /dpn{p))- For the integrable Hamiltonian H^{p) the equations of motion reduce to 
p = —dH^/dq = 0, q = dH^ /dp = \{p) with solutions p{t) = po, = <io + ^{Po)t (mod27r). Con- 

sequently the phase space is foliated by invariant n-tori with quasi-periodic motions characterized by the 
frequency X{p). The result of adding a small perturbation to the integrable Hamiltonian H^{p) is the de- 
struction of all tori except for those whose frequencies satisfy a non-resonance condition. 

In the KAM theorem attention is restricted to tori supporting quasi-periodic motions with an appropriate 
nonresonant condition. Given the set of all frequencies O C M", the frequencies satisfying this nonresonant 
condition are called diaphontine frequencies and as a subset of f2 are defined by 

I7r = {A e r2 c R"| |A- fc| > r||fc||-", VfceZ", k^Q} 

for some positive constant T . 

The method used to prove the KAM theorem in and which this paper follows, is standard in 
classical perturbation theory. The central idea is to construct a suitable canonical transformation -0 which 
brings the original Hamiltonian H into a normal form which depends on fewer angles, none if possible. 
The transformation il) is constructed iteratively as the composition of successive near-identity canonical 
transformations (p^ , (f)'^ , . . . . Taking the limit as the number of iterative steps goes to infinity results in the 
elimination of the perturbation of the Hamiltonian leaving an integrable system that is isomorphic to the 
original system. To construct this canonical transformation the authors in Q and ||^ make use of the Lie 
series method which has the advantage of avoiding any inversion and thus any reference to the implicit- 
function theorem. The Lie series method requires constructing a Hamiltonian called the generating function. 
At the fcth step of the iterative process a generating function Xk is constructed and the desired near-identity 
canonical transformation 0'^ is obtained as the flow at time 1 associated with such generating function. 
The generating function is constructed via the nondegeneracy condition and by solving a partial differential 
equation. 

The restriction to diophantine tori comes about from the problem of solving the partial differential equa- 
tion mentioned above which has the form 

j2X~iq)^Giq), (1.3) 

where the unknown function F and the known function G are defined on the torus T" and G has zero average. 

In this paper we consider two similar KAM theorems where, in each, there is an aperiodic time dependent 
term in the perturbation. In both cases we consider an n-degree of freedom real valued nearly-integrable 
Hamiltonian in action-angle variables of the form 

Hip,q,t)^H%p) + KH\p,q,t), 

where p = {pi, ...,Pn) € M", q — (qi, ...,(?„) G T", i e K and k e R is a small perturbation parameter. In 
the first case the perturbation considered has an exponentially decaying aperiodic time dependent term and 
a quasiperiodic term depending only on the angles q 

H\p,q,t) = 5fee''-«+ We'''^ 

where /(p) is a bounded function and ek{t) decays exponentially. The nature of this exponential decay is 
explained fully in the next section, (see definition 1). In the second case the perturbation consists of an 
exponentially decaying aperiodic time dependent term and a term depending in a quasiperiodic manner on 
the angles q and on time t 

H\p,q,t)= Y + E 

where f{p), ek{t) are as above and 9 = {6i, ...,9m) is the vector of basic frequencies. Although the second 
case seems more promising as far as applications is concern, we first prove the less technical first case. In the 
last section we sketch the proof of the time quasiperiodic case and state the theorem. The only difference 
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in the assumptions of the time quasiperiodic theorem involves the diophantine condition A G fir which must 
hold for the larger vector A = (A, 9) G rather than just the vector A G Z". 

In both cases we consider the nearly-integrable system generated by the Hamiltonian H{p, q, t) and prove 
a KAM type theorem where instead of tori we show the preservation of cylinders of the form T" x M. The 
proof is similar in method to the one given in §. First, the non-autonomous vector field is made autonomous 
by making time a dependent variable, t — qn+i, and grouping it with the angles. Similarly a conjugate 
variable r = Pn+i is grouped with the action variables. We write the Hamiltonian in a Kolmogorov-type 
normal form, which consists of Taylor expanding the Hamiltonian about p = 0,r = and grouping terms 
of different orders of magnitude. Note the proof of this theorem focuses on the preservation of the cylinder 
with p — 0, T — 0. This can be done without loss of generality since any other cylinder can be shifted to the 
"zero" cylinder. 

Once the Hamiltonian is in Kolmogorov normal form we seek a suitable canonical transformation which 
brings the Hamiltonian into a normal form which depends on fewer angles. The transformation is constructed 
iteratively as the product of near identity canonical transformations 0fc; i.e (f) — (f>i o (j}2 o ■ • ■ ° (j^n- The result 
is a sequence of Hamiltonians, Hi = Hq o H2 — Hi o ^2, Hn — Hn-i ° 4>n, which come closer to the 
desired normal form. Letting n ~* 00, we obtain Hoo which is an integrable Hamiltonian. The sequence of 
canonical transformations is obtained using the Lie series method. This method and other techniques from 
the theory of several complex variables will require us to extend the domain of our real valued functions and 
consider analytic functions defined on a complex domain. 

To carry out the Lie series method we introduce a generating function of the form 



Xip, T,q,t) = i- (g, t) + X{q,t) + Y^ t)P^^ 

i=i 

where ^ e M"+^. To solve for X and Yi we must solve a partial differential equation of the form 

Y.X^^{q.t)^G{q,t), (1.4) 
^ oqi 



i=l 



where Xi = ^-(0) for i = 1, ...,n, A^+i = 1, the given function G{q,t) with zero average and the unknown 
function F{q,t) are defined on T" x R. The known function G{q,t) will be an expression related to the 
perturbation and the unknown function F{q,t) will be an expression related to X(q,t) or Yi(q,t). 

For the time independent case with functions F{q) and G{q) defined on T" the partial differential equation 
can be solved in terms of Fourier coefficients Understandably though, the aperiodic nature of the time 
dependence renders the Fourier series procedure inapplicable to the problem at hand. Instead we must use 
Fourier transform methods to solve for X(q, t) and Yi[q, t) and obtain the appropriate domains of analyticity. 
Consequently F{q,t) is written as a Fourier series with time dependent Fourier coefficients fk{t) and these 
coefficients are expressed in terms of the Fourier transform of the time dependent Fourier coefficients of 



G{q,t). Solving (L4) for the unknown function F{q,t) restricts the form of the time dependent perturbation 
requiring an exponential decay in time in order to obtain convergence of the Fourier coefficients fu (t) . 

To obtain a value for ^ and completely solve for the generating function we make use of the nondegeneracy 
assumption det(dXi/dpi) 7^ i = l,...n. Note this assumption implies the frequencies are functionally 
independent and since A„+i = 1, the ratios of the frequencies A^ i = l,...,n to A„+i are functionally 
independent as well. Thus we obtain what is called an isoenergetic nondegeneracy condition. This condition 
is the key feature in solving for f . 

After applying the iterative lemma n times and taking n --> 00 we obtain an integrable Hamiltonian 
HooiP, T,q,t) in the transformed variables {p,T,q,t) which generates an autonomous vector field in T" x R 
with solutions of the form (p, f ) = (po, tq), 9 = (1 + i'iE(^A)Xt + qo, t = {1 + KE(^A)t. These solutions imply 
the phase space ]R"+i x T" X M is foliated by invariant infinite cylinders each sustaining winding motions 
identified by the frequency X{po) and evolving in the t direction. 

In Appendix ^ we present the Rossby wave flow [|l^ with a time decaying perturbation and show this 
system satisfies the hypothesis of the theorem presented in this paper. 
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2 Preliminaries 



Although we begin by considering a real Hamiltonian dynamical system, techniques from several complex 

variables will be use in the analysis. This requires complex extensions of the original functions defined in R™ 
to functions defined in C™. We begin describing these complex extensions with some notation. Denote the 
open ball of radius p centered at po by B so that po & B c M™, and we assume without loss of generality 
p < 1. Also, denote q = (qi, q„i) G where we identify functions on T"* with functions on that are 
2tt periodic in each qi, ...,qm- The complex extension of B x T™ is given by 

Dp,Po = {{P, Q) e C2™| \\p - Poll < p. Re 5 e M" mod 27r, ||Im q\\ < p}, 

where Im q = (Im qi, ...,Im and = max^ \vi\ for v € C™. Define Ap.p„ as the set of all complex 
continuous functions defined on Dp^p^, analytic in the interior of Dp^p^, 2tt periodic in the variables gi, gmi 
and real for real values of the variables. 

Consider the following non- autonomous Hamiltonian 

H{p,q,t) = H\p) + KH\p,q,t), 

where p e R", g € T" are the action and angle variables respectively, t e M and /t e M is a small perturbation 
parameter. The nature of the time dependence is explained at the end of this section. Hamilton's equations 
are given by 

. __dH _ dH^ ■ _ dH° dH^ 

^ dq dq ' ^ dp dp dp ' 

For the system generated by the integrable Hamiltonian H^{p) we define the angular frequencies \{p) = 
(^^, = (Ai, A„) and Hamilton's equations reduce to p = —dH^/dq = 0, q = dH^ /dp = A with 

solutions p{t) = po, q{t) = qo + \{pQ)t. Consequently the phase space is foliated by invariant n-tori, each 
of which sustains quasi-periodic motions characterized by the frequency A(po)- 

To analyze the full Hamiltonian H{p,q,t) with the time dependent perturbation H^{p,q,t) we change 
this non-autonomous system into an autonomous system by making time a dependent variable and effectively 
extending the phase space from R" x T" to R" x R x T" x R. For notation sake we group t with the angles 
q and its conjugate variable r with the actions p and write the Hamiltonian H{p, r, q, t) with {p, r, q, t) = 
y,q') € R" X R x T" X R 

H{p', q') = H{p, T, q, t) = H\p) + kH\p, q,t) + T 

or 

H{p', q') = H'^ip, r) + H\p, q, t) = H\p') + kH\p' , q'), (2.2) 

where II°{p') = H^{p) + t. Given 0<a<p<l,wc define the complex extension of M" x R x T" x R 

Dp,a,po = {{p,T,q,t) e C2"+2| \\p-po\\ < p, \t\ <a, Re g € R" mod 2w, ||Im q\\ < p, jlm t\ < a}, 

and define Ap,a.pa as the set of all complex continuous functions defined on Dp^^^p^, analytic in the interior 
of Dp^cr,poj 27r periodic in the variables qi, ...qn, and real for real values of the variables. Hamilton's equations 
for the new autonomous Hamiltonian are given by 



— -M. — .- dH^ 

(2.3) 



~ dt ~ at ' 



^ ~ an ~ an ~^ an ■ 



dH _ -I 
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Of the new Hamilton's equations, ( |2.3| ), the first, third and fourth equations are equivalent to the original 
system (2.1). For the system generated by the Hamiltonian H^{p') Hamilton's equations reduce to p' = 
-dH"/dq' = 0, g' = dH°/dp' = X{p') with solutions p'{t) = p^, q'{t) = A(p(,)t + 50 where 

We will often write A = (A, 1) = (Ai, A„, 1). It is instructive to write out the solution q'{t) — 
{qi{t), ...qn{t),t) — (Ai(p'), A„(p'), l)t + qo where go G x K and the 71 + 1 term of qo is zero. Con- 
sequently the phase space {p',q') G (R""'"^,T" x R) is foliated by invariant infinite cylinders each sustaining 
winding motions characterized by the frequency X{po) which evolve along the t direction. 

For a vector v G C" we define the norm \v\ = \vi\. For a scalar valued function / e Ap^a,po we use 
the norm, ||/||p,o-jpQ = sup^^/ j^/^gj^^ \f{p',q')\. For vector valued functions / — (/i, f2n+2) with values 



-■2n+2 define / G Ap^a,po if fi G •A.p,a,po^ [i = 1) 2n + 2), and the norm is defined as follows 



= max II/, 



\\J llp.O-.Po — ^^^'^^ ||./j||p,<T,po- 

% 

For a (2n + 2) x (2n + 2) matrix C whose entries C^j belong to Ap^cr.po we write the norm of C as ||C||p_CT,po = 
maxjj ||Cij ||p,cr,po- -'-f t^*^ matrix C is a constant matrix then ||C|| = max.ij- |Cij |. 
For the integrable Hamiltonian recall the frequency map 

The set of all possible frequencies in R" is denoted by 17 = \{B). The subset of diophantine frequencies which 
will be important throughout the paper is defined for some positive constant F by 

nT^iXeVtd M"||A • k\ > r||fc||-", Vfc e Z",fc ^ 0}. 

Given a function f{p, r, q, t) = /(p', q') e Ap^a periodic in each of the qiS except qn+i — t, f{p'*) will 
denote the average of the function / over the angles gi, g„ and evaluated at p'* 

7(p'*'t)^7^ r ■■■ r f{p'\q:t)dq^---dqn. 
\.^^) Jo Jo 

Given a function f{p,T,q,t) ~ f{p',q') G Ap,a 27r periodic in each of the g^'s except qn+i = t, we will 
denote the average of the function f{p', q') over the angles qi, qn and over t evaluated at p'* as follows 

— 1 / 1 ^ " '"^'^ ''^'^ '""^ 



2TUj X ■■■Jo (2-4) 

Note for a given Hamiltonian H e -4p.o-,po: one can assume ||i/||p,cr < 1- This can be obtained by the 
change of variables {p',q') — > {ap',q') with a suitable positive constant a. 

Now we address the topic of the time dependence of the perturbation H^{p,q,t). Since we want to 
consider a time aperiodic dependence and the perturbation plays a key role in the solution of the canonical 
transformation, the perturbation has to have a particular form for certain Fourier transforms to converge. 
The form of the perturbation consists of an exponential decay in time of the time dependent part of the 
perturbation. The exponential decay in time of the perturbation is described by the following definition. 

Definition 2.1 



We say a complex valued function f{z), z = x + iy with x,y G R, is of {Ci,C2,ci,C2,h', IJ.)- exponential 
order with respect to x and write 



O {er^^'^^y-') [x -> 00) 
O (e(^-^'^) {x -00) 
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where v,iJL,e>Q,v>e,ii>e if there exist constants Ci, C2, ci > 0, C2 < such that 

\f{z)\ < Cie-f'^-^)^ for < ci < a; < 00 and 

< Cse^''-^)^ for- 00 <x <C2<0. 
The iterative lemma will require the perturbation to have a particular form. This form is defined below. 

Definition 2.2 

We will refer to a function F{p',q') as having {Ci, C2, ci, C2, v, iJ,)p' , q' -exponential form if it can be ex- 
pressed as follows 

F{p',q')=f{p')+r{p',q')+g{p',q'), 
r{p',q')= E 

9{p',q')= E hk{p')ek{t)e^''-^, 



where F{p' , q') G Ap,a- and ek{t), t — tj^ + itj, tji, tj G M, is of Ci, C2; ci, C2- exponential order with respect to 
tji. When it is not important for context to specify the constants Ci,C2,Ci,C2,v, n we will simply refer to 
functions of exponential order with respect to x or functions of p' ,q' -exponential form. 

3 Statement of the Theorem 

We now state the main theorem in terms of the Hamiltonian H{p' , q') = H^{p) + kH^{p, q, t) + r. 
Tlieorem 3.1 

Consider the Hamiltonian H{p',q') = H°{p) + KH^{p,q,t) + r = H°{p') + KH^{p',q') of 
{Ci,C2,ci,C2,y,n)p' ,q' -exponential form defined on B x where H^{p') = H^{p) +t . Fix p'q G B and 
denote 

Assume there exists positive numbers F, 7, p, a, d, k, all less than one and L, such that p> a and 

Aef^r, |A|<i, 



.PO' 

-lIL-.l 



d\\v\\ <\\C*v\\ <d-^\m, VweC", 
moreover, ||p,CT,po < 1> T^^^ last condition will imply the isoenergetic nondegeneracy condition on I). 

Then there exists positive numbers E, k' , E' X' , f,C3,^, p' and a', such that if ||-f?^||p,(7,po — ^ one can con- 
struct a canonical analytic change of variables 

ijj '. D pf ^fj' ^ Dp^fj^pQ^ 

{V',Q') ^ i'{V',Q') = {p',q'), 
with tp e Ap'^a', which brings the Hamiltonian H into the form H' = H o tp given by 

H'{V', Q') = {H o i;){r', Q') = a' + (1 + k'E'C)X ■ V + ), (3.1) 
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where a' G K. The change of variables is near identity in the sense that 

II - identity\\p>^a ^ as H-ff^llp.o-.po ^ 0- 

In particular, we can take 



In the new coordinates Hamilton's Equations are given by 

r) H' 

Q' = g^pjiV, Q') = (1 + n'E'OX + 0{\\V'\\), 
V' = -^{V',Q') = 0{\\Vr)- 

The solutions are given by 

Q'{t) = (1 + n'E'OXt + Q'o V'{t) = P^ = {Vo, fo) = 0. 

We can write out explicitly these solutions Q'{t) = (Qi(t), Q„(t),T(t)) = (1 + n'E'O^t +&„ = {! + 
k'E'('){X, l)t + Q'q which indicates in the new coordinate system the phase space R" x R x T" x Ris foliated 
by invariant infinite cylinders each sustaining quasi-periodic motions identified by the frequency X{Po) o-nd 
evolving in the t direction. Note under the transformation the time variable T{t) = (1 + k' E'C,')t has shifted 
by a constant proportional to the size of the perturbation. 

4 The Normal Form 

We begin the proof with a trivial rearrangement of the Hamiltonian that consists of Taylor expanding in p, r 
around p' = 0. Recall the Hamiltonian 

H(p',q') = H\p)+T + KH\p,q^t) = H\p') + KH\p',q'). 

Taylor expanding the Hamiltonian H{p' ^ q') about = we obtain 



H{p',q') = if(0,90 + |^(0,9')-p'+2|^(0'«')(p',p')+O(lb'll') 
= F°(0) + kH^{Q, q') + \-p' + Kdp,H^{Q, q') ■ p' 

+ \dl,H\0){p',p') + K\dl,H\i), q')ip',p') + 0{\\p'\n 
The terms can be grouped by order of magnitude in p' 

, n+l 

H{p', q')=a + X-p' + A{q') + B{q') • + 3 E ^i^P'iP'o + q'), (4-1) 

where 

a = H(0) = H0{0) + kW{0) = H°{0) + kW{0), 
KA{q') = H{0,q')-a = K{H\0,q')-W{0)), 
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and a G M, A, Qj, i? € Apu,R = 0{\\p \\). This Hamiltonian has the desired form (3.1) of Theorem 



except for the terms A{q') and B{q') which constitute the actual perturbation of the Hamiltonian ( [l.lD 



5 Bounds on the Normal Form in terms of the Original Hamilto- 



nian 



I n w hat follows we present several results that allow us to reformulate Theorem |3.l| for the Hamiltonian 
( [4.l| ) in normal form. Since we are interested in estimates for A,Bi,Ci_j and these are defined in terms of 
derivatives of the original Hamiltonian, we often make use of Cauchy's inequality. For scalar valued functions 
/ G Ap^cr, S < a < p and nonnegative integers fc^, Z^, i = 1, n + 1 



^fci+...fe„+i+ii+...+i„+i 



fip,q) 



< 



fci!...fc„+i!/i!...l„-|-i! 

jfei+...+fc„+i+Zi+. ..+;„+! 



\f{p, q) e Vp^sM-S 



Lemma 5.1 

Recall the Hamiltonian of interest H(p',q') = H^{p) + kH^(p, q,t) + t = H'^{p') + nH^ijp' ,q') where 
H^{p') = H^{p) + T. Given the n x n matrix 

assume there is a positive number d < 1 such that d\\v\\ < \\C*v\\ < d^^\\v\\ \fv G C". Define the (n + 1) x 
(n + 1) matrix 

it follows ||C*f|| < Vu G C"^"'^. Proof See appendix [B|. 



Lemma 5.2 (Lemma on bounds) 



Given a constant E such that \\H^\\p.a < E 



1. max{\\A\\p,,,\\B\\p^,) < 2f = Ei. 



2. There exists a positive number m < 1 such that \\Cv\\p^a- < "m ^\\v\\, \fv G C"^"'^. In particular we 
take m — ^. Proof See appendix 

Lemma 5.3 (Isoenergetic Nondegeneracy Gondition) 



Assume Y.\K\ < L. Given \\C*v\\p^a > rfNI, V5 G C", 2m > L. It follows \\lv\\p^a > l\\v\\, Vu G C"+2 
for some I = ^\d — L\. Proof See appendix 
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Lemma 5.4 

Given ||Iu||p^cr > ^Utillpjcr with I — ^\d — L\. It follows 

A 



>f\Hp,a, 



where / = |Z- 
Proof See appendix 

6 Re-statement of the Theorem for the Normal Form 

One can deduce, based on the bounds on the normal form in terms of the original Hamiltonian, theorem | 
from the following theorem. 

Theorem 6.1 Main Theorem 

For given positive numbers K,T,"f, p,a,m all less than one and L consider the Hamiltonian H{p',q') 
U{p' ,q') + P{p' ,q') of {Ci,C2,Ci, C2 J V, n)p' , q' -exponential form defined in Dp ,j by 



P{p',q') = n(^Aiq') + J2B^{q')p'^ , 



with \\H\\p,„ < 1, where A G fir and |A| < L,A,Bi,Ci_j,R G Ap^a o,nd R is of order . Assume 

2to||w|| < 11(7*^11, VueC", 

\\Cv\\p^„ <m-^\\v\\, VweC"+^ 

Then there exists positive numbers Ei, Koo, Ej^ , Coc, p' , jin' , f all less than one and L',C3,j also positive 
with 

p < p,a' < (7, m' < m, 



f\\v\\ < 



C X' 



such that if 

max(||A||p,„,||B||p,^) < El, 
we can construct a canonical analytic change of variables 

with (f) G Ap'^a' , which brings the Hamiltonian H into the form 

H\V, Q) = {Ho <j>){V, Q) = a' + X'-V + R{V, Q), 
where H' is of p' , q' -exponential form, A' = (1 + Koo-Ef°Coo)A, a' G M and R ^ 0(||P|p) G Ap'. 

The change of variables is near the identity in the sense that 

||(/) — identityWp' ,a' as \\P\\p^a 
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In particular, we can take 



In the new coordinates Hamilton's Equations are given by 

r)H' 

Q' = ^{V, Q') = (1 + Kooi^rCoo)A + 0{\\V'\\), 

r' = -gQ^{v',Q') = o{\\v'\n 

The solutions are given by 

Q'{t) = (1 + /«oo^rCoo)At + Q'o V'{t) =V'o = {Vo, to) = 0. 

We can write out explicitly these solutions Q'{t) = Q„(t),T(t)) = (1 + KooEi°Coo)^t + Q'q = 

(1 + KooEi^C<x:)i^, i)t+ Q'q which indicates in the new coordinate system the phase space R" x R x T" x M is 
foliated by invariant infinite cylinders each sustaining quasi-periodic motions identified by the frequency \{Vq) 
and evolving in the t direction. Note under the transformation, the time variable T{t) = (1 + Koc-Ei^Coo)i 
has shifted by a constant proportional to the size of the perturbation. 

7 Perturbations tending to g— Periodicity and t Quasiperiodicity 

We arc considering two similar KAM theorems each involving an n-degree of freedom real valued nearly- 
integrable Hamiltonian in action-angle variables of the form 

H{p,q,t) = H\p) + KH\p,q,t), 

wherep = (pi, ■■■,Pn) G R", q = {qi, g„) G T", t eM. and k is small. After extending H to the complexified 
domain Dp^cr we have H G Ap^a- In the first case the perturbation considered has an exponentially decaying 
aperiodic time dependent term and a quasiperiodic term depending only on the angles q 

H\p,q,t) = 9ke'''+ E fkip)ek{t)e^''-^ 
feez" feez" 

where fk{p) is a bounded function and ek{t),t = ta + itj, is of exponential order with respect to t^. In the 

second case the perturbation consists of an exponentially decaying aperiodic time dependent term and a term 
depending in a quasiperiodic manner on the angles q and on time t 

H\p,q,t)= J2 + ^ A(p)efe(t)e*«, 

where fk{p), ek{t) are as above and = {6i, ...,0m) the vector of basic frequencies. 

For each of the cases above, it is necessary for the application of the Lie series method to construct a 
generating function. For this purpose we will formulate an expression in terms of the perturbation H^{p, q, t) 
evaluated at j> = , call it G{q,t), and solve the partial differential equation 

"+1 op 

Y^X~{q,t) = G{q,t). (7.1) 

2—1 

For the first case, the given function has the form 



G{q,t) = GA{q,t)+GQ{q) = ^ fffe(t)e^'=-« + ^ ~gie''-' 



feez" 
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and we look for a solution of the form 



Fiq,t) = FAiq,t) + FQ{q) = ^ /fc^e*' + ^ 



il-q 



Substituting in (7.1) we obtain 



dF 



n+l 

i=i ^* fceZ" 



i\-kfk{t) 



Clearly — G = becomes 



E 

( [7.;^ is satisfied if 



dfkjt) 
dt 



Jk-q 



li-q 



(7.2) 



*(A-fc)/fc(t) + %^=5fe(t) 

i(A-/) 

and the problem reduces to two problems. Similarly for the second case, the given function has the form 
G(g,t)-G^(g,t) + GQ(g,t)= ^ .g,(t)e'^-'+ ^ 

we look for a function of the form 

F{q,t)^FA{q,t) + FQ{q,t)=Y,Mt)e^''-''+ ^ fi^'<''^''\ 



feez" 



and the problems to be solved are 



fi 



' dt 

9l 

i(\,0)-l 



Therefore in both KAM theorems solving ( [m| ) is equivalent to solving two p.d.e's; one p.d.e involving 
quasiperiodic dependent functions, (GQ{q),FQ{q) or GQ{q, t), FQ{q, t)) and a p.d.e involving p', g'-exponential 
form functions ( GA{q,t), FA{q,t)). In the next section we give the result for solving the p.d.e involving 
quasiperiodic functions. 



8 P.D.E on a Torus 

Consider the real valued functions F{q) and G{q) where q = {qi, qn) £ C" and Re g = (Re qi, Re g„) € 
M" mod27r. That is, F{q) and G{q) are functions on C", 2tt periodic with respect to the real part of each 
gi,...,(7„. Although we are concerned with real Hamiltonian systems, we will use techniques from several 
complex variables theory in the analysis. This will require complex extensions of functions originally defined 
on R" to C". Fix p < 1, we define the complex extension of T" 

Dp = {qeU'\ Re g e K" mod 27r, ||Im q|| < p}, 

where Re q = (Re gi, Re Im q = (Im gi, ...,Im g„). We define Ap, the set of all complex, continuous 
functions defined on Dp that are analytic in the interior of Dp and real for real values of the variables. 

The following lemma gives a useful bound for the Fourier coefficients of a function G(g) with a 2tt periodic 
dependence on g. 
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Lemma 8.1 



Assume G e Ap and G{q) = Z^fcgZ" fffee*'''^, then for every G Z" \gk\ < \\G\\pe ^''^p, where \k\ = 
\ki\. Proof see appendix 

The following lemma gives the existence and uniqueness of solutions of a p.d.e involving quasi-periodic 
functions. The lemma also gives analyticity results for the solution, F{q), of the p.d.e given analyticity 
restrictions on the known function G{q). 

Lemma 8.2 Consider the following linear partial differential equation 

Y^X., — {q) = G{q), (8.1) 
where F and G are functions defined on the torus T", and assume A = (Ai, A„) G Qr, for some F > 0, and 



G G Ap for some positive p < 1 with G = 0. Then, for some positive 5 < p, (8.1) admits a unique solution 
F G Ap-s with F = 0, and one has the estimates 



||F||p-.<^||G||p, 



dF 



dq 



p-d 



where vj = 2''"+^ P''~oof see appendix 



9 P.D.E on a Torus with Time Dependent Coefficients 

Consider the functions F{q,t) and G{q,t) where q = (qi,...,g„) G T", i G M. We will require as before 
complex extensions of functions originally defined on R"+^ to C""*"^. We define the complex extension of 

T" X R 

Dp,a = {iq,t) G C"+i| Re g G M" mod 27r, ||Im q\\ < p, |lm t| < a}, 

where Re q = (Re qi, ...,Re g„), Im q = (Im qi, ...,1m </„). We define Ap^a, the set of all complex, contin- 
uous functions defined on -Dp, a that are analytic in the interior of Dp ,j and real for real values of the variables. 

We want show that for a given function G{q, t) G Ap^a satisfying certain conditions there exists an F(q, t) 
which satisfies the following 

Y.\^{q.t)^G{q,t), (9.1) 

where we identify qn+i with t and set A„+i = 1. 
Theorem 9.1 

Given 



Y,\~{q,t)^G{q,t), 
^ oqi 



Giq,t) = J2 5fe(^)e^'■^ F{q,t) - ^ f^it) 



where G G Ap,a for 1 > p > 0, 1 > cr > and q G T",t = tji + itj G C,tn,ti G R. Since G G Ap.a it 
follows gk{t) is analytic in the strip < a and further assume for any given k 



9k{t) 
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where z/, /i > 0. Then there exist S > and ^>Q,^<v,^<pL such that fk{t) is analytic in the strip 
\ti\ < a — S and satisfies 



h(t) = 



0(p(M-£)tH)c)(e-|fclp) (t^ ^ _oo) 



Proof 

Since F[q^t) and G{q,t) are 27r periodic in each qi, ...,(7„, one can write the Fourier series 



F{q,t) = J2 fk{t)e^^-\ Giq,t) = ^ g^^e^'^-' 



Substituting these in (9J) results in the following differential equation which the Fourier coefficients fkit) 
and gk{t) must satisfy 

iCX-k)Mt) + ^^gk{t). 
The solution of the differential equation is given by 

hit) = /fc(0)e-*(^-^-)* + /* gkis)e^^'-'^^^-'Us. (9.2) 
Jo 

We can argue fk{t) has a Fourier transform provided /fc(0) satisfies certain condition. Ultimately, with this 
condition on /fe(0), we will obtain a new expression for fk{t) which will be better suited to apply Fourier 
theory. Assuming that gk{t) has a complex Fo urier transform ^(cj), with u; = u + iv, u,f G M, analytic in 
some strip — < — /3 < v < (3 < v, rewrite (|9.2|) in terms of the Fourier inversion formula. That is, given the 
Fourier transform 

gk{t)e~'^'dt, 



Gk (t^) 



(9.3) 



and the Fourier inversion formula 



9k{t) 



oo+i/3 



gk{io)e'^'duj, 



becomes 



/fe(<) = /fe(0)e-'(^-'=)* 



oQ-\-i(3 



oo+i/3 



gk{uj)e''^'duj e*(^-'=)(*-*)ds. 



Interchanging integrals and integrating with respect to s gives the expression 



/fc(t) = /fc(0)e-^(^-'='* + e-«(^-'^)* 



Qk{uj) 



pi(u>+\-k)t _ 2 

i{uj + A • fc) 



(9.4) 



We consider an arbitrary function s{lj) with the property J^^^^p s{uj)dLu = 1 and rewrite (9.4) as follows 



fkit) = e-'^^-'^^* 



oc+i/3 



fkiO)siuj) + GkiLo) 



oo+'i/9 



— oo+'i/3 



/fc(0)5(c.)e-'(^-'=+")* + gfc(L.) 



gi((.j+A-fc)t _ 2^ 

i{uj + X ■ k) 
i{uj + X ■ k) 



dio 



(9.5) 
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Note ( p.5[ ) is written in the form of the Fourier inversion formula. That is 

(>QO-\-iP 



fkit) 



00+1/3 



(9.6) 



where J-'k the Fourier transform of fkit) given by 



Consequently, from (9.5) and (9.6) we write 

M^) = /fc(0)s(^)e-^(^-'=+")* +5fc(a;) 



2 _ g-i(ii'+A-fc)t 



i{uj + X ■ k) 

Furthermore, note J^k is only a function of lu which implies the following 



(9.7) 



/fc(0)s(c.)e-'(^-'=+")* = gfc(c.)- 



, — i{uj+X-k)t 



i{uj + X ■ k) 



or 



fk{0)s{Lu) 



Gki^) 



i{LU + X ■ k) 



(9.8) 



(9.9) 



and (9^) reduces to the following 



'-oo+i/3 i{uj + A • k) 

Note the condition on /fe(0) given by (9J) reduces to ( j.lC| ) for < = by integrating (9^) on both sides 



(9.10) 



/fc(0)s(tj)dcj = / TTT'^'^' 

-oo+'i/3 J-oo+i,3 «(W + A • fc) 



^fc(c^) 



-oo+ip J -oo+ip lyuj + X ■ K) 

/,(o) . g^'y' 

J-oo+ip i{uj + X - k) 



Now, using Lemma B.3 it follows 



0(g-(i.-e)t«)0(g-|/c|p) (i^^oo) 



5fe(i) = 



or 



0(e(M-£)tH)C)(e-|fc|p) (ifl ^ -oo) 
\9kit)\ < Cie^(''-^)*«e-l'=l'', < ci < < oo, 



l5fc(0l < C2e(^-^)*«e-l'=l^ -(X3 < < C2 < 0. 
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where the k order will pass through the proofs of Theorem E.l and Theorem E.2 like a constant. By Theorem 



E.l, there exists a, S > such that the complex Fourier transform of gk{t) defined by 



is analytic in the strip —fi^6<v<h' — S and satisfies 



CC3 



CC3 + 



iS-e) 



where we have used the fact 



iS-s) 



, < u < oo, 



-oo < u < 



and C3 = min(ci, |c2|), C3 — max(Ci, C2). From (9.7) and ( p.8\) we see that the complex Fourier transform 
of fkit) is given by 



i{X ■ k + Lj) 

We will show that J-k{uj) is exponentially small with respect to u and k and analytic in two parallel strips. 
Clearly, since A • A: e R, we can pick a 7 > such that |A • A; + w| > 7 and since Gkii^) is analytic in the strip 
—fi < V < v, J^k{^) will be analytic in the strips ^(7 < v < v) and B{—^i < v < —7). We have 



fkit) = 



oo+z/3 f \ 



-00+1/3 i{X - k + uj) 

and the integral converges uniformly for \tj\ < <j. Hence fk{t) is analytic in this strip. Furthermore 



fk{t) 



and 



\fk{t)\ < 



-00+1/3 i{X - k + u;) 
\gkiu + iP)\ 



-00 i(A • fc + w + i/3) 



00 ■ k + u + il3)\ 



i{u+tl3){tR+iti) 



du 



< - 



< 



\gkiu + iP)\\e 



iutRll ~uti I |„-/3tH| I -i/3t 



27r7 
2 



2n'-f 



f,-\k\p^-l3tR 



Cc. 



Cc. 



{6-8) 

C3e-(«-e)« 



(6-8) 



— (o- — e + </) (o- — e — </) 



> 



where we have used the fact that we can pick a 7 > such that • fc + u + i(3)\ — (A • /c + u) + /3 
/3 > 7. For the u intervals [0, 00) and (—00, 0], we pick tj = —a+5 and t/ — a~5, 5 > 8 , respectively, so that 



1/2 



l/fcWI < 



4 e-^^^Pe 



k\p -0tR 



27r7 {S - 8) 



C3C3 + 



iS-e) 
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Finally, we obtain the order results by taking /?, which was arbitrary from the definition of G{i-u), arbitrarily 
close to v and —fi. Note that since fk{t) is only defined in the strips ^(7 < v < i^) and /i < v < —7), /3 
must be within these strips; ^ < (3 < v and — /i < —(3 < —7. By Lemma A.2 for some 6 < p 



\F\ 



P-^^^-^ " 7\/2^ (<5 - e) 



(<5-^) 



We can easily find estimates for the partial derivatives of F{q,t). Given 

F{q,t)= J2 fk{t)e''"''^ Fg^p,,, 



the partial derivatives of F{q, t) with respect to qj j ^ n + 1 are as follow 



dqj 



For some positive S < p 



1 



\tk,Mt)\ = |fc^.||/,(t)| < ( 4)elfclV,(t)| 



< 



,\k\5 



4 e-l'^l'' 



CC3 



(J -(S-e)cs 



{6-6) 



By Lemma A.2 with p — (5 in placed of p and (5 in place of S we obtain the following 



dF 



< 



4 1 1 



CC3 



and e ^ 



Similarly, we obtain an estimate for the partial derivative of F{q,t) with respect to 



9qi ~" " >-i5-(5,CT 

qn+1 — t which is given by the following expression 



^ dfkifl ik.q 
dt ^ dt 



Using Cauchy's Inequalities 
dfk{t) 



dt 



1 1 4 e-l'^IP 

- 5ll-^'=ll^-^7V2^(<5-e) 



CC3 + 



iS-s) 



By Lemma A.2 for some 5 < p 
dF 



dt 



1 4 1 
< -- 



p-5,a-2S 



6 7\/27r [S - e) 



CC3 + 



(<5-£) 



and 



(9F 

for « = l,...,n+L 



Let d = max(l/e(5, 1/(5), then 



dq' 



< 



Ad 1 



C7c3 



□ 
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10 Iterative Lemma 



Lemma 10.1 (Iterative Lemma) 

Given positive constants T,p,a,K < 1, consider the Hamiltonian H{p',q') = U{p',q') + P{p',q') defined 



by 



U{p', q') = a + A • y + i ^ {q')p[p'^ + R{p' , q'), 
Pip',Q') =>^^A{q')+^B,iq')p'^, 



with ll-ff llp^o- < 1, -ff G Ap,a, with some constant vector A — (A, 1) G C"+^ such that |A| < Lq for some constant 
Lq, a e Or, A, Bi, Cij,R G Ap^a end R is of order \\p\\^. Assume H is of (Ci, C2, ci, C2, y, n)p', q' -exponential 
form. Assume there are positive constants m, f,Ei<l such that 

2m||{;|| < lie*!;!!, VS e C", 



fM < 



C A^ 



A 

\\Cv\\p,,<m-^\\vl Vi;eC"+\ 

max{\\A\\p,a, \\B\\p^c,) < E\. 

One can construct positive constants pt,,at,,mi,, L,C3,j with < p, cr* < a, m* < m, L = |Lo;C3 
niin(ci, |c2|) and 

A 



V = 



kEi, K = 10n + 9, A = 



3"2^»+iW(6n + 6)4(c3 + l)\ .+, _ .^n+i fn + l^-+' 



such that for any 6 > small enough so p ~ AS > p^,, cr — 4^ > and with kEi small enough that 



m ■ 



4(n + l)r] 



>m*, / 



(n + 3)77 



(10.1) 



there exists a G M and an analytic canonical change of variables, : Dp_4s,a-4S ^ Dp.a, 'P G Ap-j^s^rr-AS, 
which transforms the Hamiltonian H to H'{V', Q') =UH = {H o (j>){V' , Q'). H'{V', Q') is of V, Q' ' expo- 
nential form and can be decomposed as the original H with primed quantities A',B',C', and R' 

H' = U' + P', 

i,3 



where 



P' = k' [A'iQ') + B'{Q') ■ V) , 



A' = (1 + k£iCa)A G Qr, a' = H'{0), A'{Q') = H'{0, Q') - a', 

f)H' rP'H' 

BKS') = ^(0,S')-A:, C[^{Q') = ^^{Q,Q'l R' = 0{\\V'\\% 
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and satisfies similar conditions with positive numbers p' , a' , m', /', k', E[ all less thans one and L' also positive 
with \\H'\\p, < 1, where A' , B'i,C^j, R' G Ap'^a', 



f'\M < 



C X' 



-<n+l 



where 



A 

\\C'v\\p>,^, <m'-^\\v\\, Vt;eC"+\ 
max(||A'||p',^', ||-B'||p',a') < E[, 



I as: / as: > 4(n + , (n + 3)?? ^, 3^ 
p'=p-4S>p^, C7 =(T-4(5>(T*, m' = m — > m*, f = f-^ L' = -Lq, 



at 



at 



A 



:k\ E[ 



A El 



Furthermore, one has for any function F € Ap^a, — F\\p',a-' ^ v\\P\\p,o-- 
Proof 

Wc will construct a canonical transformation of the Hamiltonian, with a generating function for this 
transformation denoted by x- The generating function is constructed in such a way it eliminates the part 
of the transformed Hamiltonian that prevents the preservation of the invariant structure p' = 0. Recall the 
Hamiltonian 



H = U + P, 

Uip', q')=a + X-p' + lYl CidWiV'j + Rip', q'), 



P{p',q')=^(^A{q')+Y,Bi{q')p'^ . 



We write the transformed Hamiltonian as follows 

H'=UH = U + P + {x,U} + [{x, P}+UH-H- {x, H}] . 
Alternatively, assuming analyticity in t and taking the Taylor expansion, one has the Lie series 



UH=y^ —.L^H, 



m=0 



where we denote i^-ff = Fl and L'^H = {L^ ^H,x} for m > 1. For the wth remainder of the Lie series, 
we use the notation 



1=0 l=m 

With this notation, the transformed Hamiltonian can be expressed as follows 

H' =UH = U + P + {x,U} + r^iU, x, 1) + ri(P, x, 1), 
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so that {x, P}+UH-H- {x, H} = r-iiU, x, l) + ri{P, x, !)■ We choose x so that {x, P}+UH-H- {x, H} 
is second order in the size of the perturbation and 



P + {x.U} = c\-p' + 0{\\pT) 



(10.2) 



which iniphcs P + {x, U} does not contribute to P' but it shifts the frequency by a multiple of the original 
frequency. Furthermore, one can show all conditions are met by a generating function, similar to the one 
introduced by Kolmogorov, of the form 



X^X{q')+S,-q' + Y{q')-p\ 



(10.3) 



where the functions X{q') and Y{q') are of g'-exponential form and £^ G M"+-'^. A simple calculation with the 
appropriate definitions of U and P gives 



{X,U} = 



Equivalently 



Y.\i^{q')p[ + o{y\ 



(a • C + A • d,,X{q')) + ( - e • Ciq') - d,.X{q') ■ C{q') - A • d,,Y{q')) ■ p' 

+ OiWp'f)- (10.4) 



To obtain the form of (10.2) we first impose the following condition 

- \ ■ C - \ ■ dq,X{q') + KA{q') = 0. 



We set 

so that dlOSl) becomes 



A • ^ = kA, 
X-dg,X{q') = K(^A{q')-Tj 



(10.5) 
(10.6) 

(10.7) 



By Lemma F\2 and the definition of A{q') we know the right hand side of ( 10.7 ) consists of a quasiperiodic 
part and a part of exponential order with respect to time. Consequently ( 10.7 ) can be split in two problems 
corresponding to the quasiperiodic part and to the part of exponential order with respect to time respectively 
and solved as ex plained in section |^. By lemma 8^ and theorem |9.l| there exists an analytic function X{q') 
satisfying (10.7). Note, by assuming the perturbation of the Haniiltonian, mainly H^(p',q'), consists of a 
quasiperiodic part and an exponential-order- with-respect-to-time part, we are assured X and dqiX are of 

the same form. Next, for some {Ca,£,) G K x R"+^ to be determined, set {kB — C ■ ^ — C ■ dq'X) = Ki?iC^A. 

Consequently we have the following conditions C • ^ + k-EiCaA = kB — C ■ dq>X, X ■ £, — kA, or equivalently 



C A^ 
A 



kB-C ■ dq'X 



(10.8) 



It follows from lemmas 5.3 and 5.4 there exists a solution, (^, C^a) G M" x M, for the matrix equation above 
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Up to this point we have determined conditions on X{q') and (^a) G M""*"^ x R to obtain the Hamiltonian 
H'ip', q') = U + kCaX ■ p' + OiWp'f) + UnBiq') - C{q') ■ ^ - C{q') ■ d,,X{q')) 



where we have set 



~{kB~C -^-C ■d^)\ -p' +^-\-d,,Y{q')yp' 

+ [{x.P}+UH-H-{x.H}] 

U + P-p' +(-\- dq,Y{q')] ■ p' + Ua, 



U = U + kE,Ca\-p' + Oi\\p'\\^) 

= a + (1 + kE,Ca)X -p' + IT. C^MM + ^(P'' l') + 0(lb'll') 



^,3 



where is 0(b'lr),/3= {nB{q')-C{q')-i-C{q')-dq,X{q'))-{KB-C-i~C-dq>X) , A' = {I+kE^Ca)^ 

and 

nA^{x,P}+UH-H-{x,H}. (10.9) 



Next we set — \- dq'Y{q') + • p' = or equivalently 

A • d,,Y{q') = p. 



(10.10) 



Note, by the definition of B and C and by the form of dg'X, ( 3 consists of a quasiperiodic part and a part of 
exponential order with respect to time. Consequently, (lO.lC) can be split into two problems and solved as 
explained in section |^. The final Hamiltonian after one application of the generating function is therefore 



H' =UH = 



^'■p'+oT. c^A<i')M + WW')] +nA = a + \'-p' + o(b'||2) + Ua. (10.11) 



where TZa is 0{k^). By lemma F.8 the new Hamiltonian H'{p', q') is of p' , (/'exponential form. The formula 
for X{q') can be obtained in terms of Fourier coefficients. Consider the following 



H\p' ,q,t) ^ G{p' ,q) +T{p' ,q') e Ap, 



(10.12) 



G{p',q) = 4(p')e^'-^ e Ap, T{p',q') = hW)el(t)e^''-' e 



X{q,t)^y{q)+T{q'), y{q) = ^ yfee^'=■^ r{q') = J] Xk{t)e^'^-'^ . 



We have assumed e\.{t) is analytic in the strip —a < Im (<) < a and satisfies 

C)(e(M-£)tfi) (tR -oo) 



e^i) 
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or 



for some positive Cq, Cq and 



< C^oe^^""^*^, -oo < t^, < C2 < 0, 



\hl{0)el{t)\ < Cie-('^-")*«e-l''^''', < ci < tj? < oo, 



\hl{0)elit)\ < C2e('^-^)*«e-l'=^l'', -oo < tR < C2 < 0. 
Recall the equation to be solved 

X-dg,X{q')^K(^Aiq')~T] 

and by the defined normal form 



(10.13) 



nA{q') 



, fc6Z"\0 



Note A = Q. Solving ( 10.13 ) can be done by separating the quasiperiodic part and the exponential-time- 
dependent part and solving the two parts separately as was presented in section four. The quasiperiodic part 
is as follows 



A • d,,y{q) = K (g(o, q) - G(0)) = ^4(0) 



Since we assume A e fir and since the right hand side of ( lO.l^j ) has zero average, by lemma 3.2 
positive 6 < p , (10.14) admits a unique solution 



(10.14) 



for some 



y{q)= Y y^^e*' G -4p_5, 



with 3^ = and one has the estimates 

k-|<^(^)"||(G-^)(0)|| 



'\k\{p-S) 



(10.15) 



dy 



dq 



p-5 



<^^\\iG-Gm\\p, 



\\y\\p-s<Y6^\\iG-Gm\\p, 

where uj = 2'*"+^ ^^'^ exponential-order- with- respect-to-time part of ( 10.13 ) is as follows 



\-d,,Tiq')^KTiO,q'), 
which reduces to the differential equation 



(10.16) 



By Theorem EA and lemma EA , there exists a (5 > such that the complex Fourier transform of e^(t) 
defined by 

-oo+ip 
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exists and is analytic in the strip —fi -\-S<v<iy — S and the Fourier inversion formula is 

/oo 
-oc 

We found, ( ^.ICl ), Xk{t) has the following form 



By Theorem 9.1 there exist (5 > 0, (5 > 0, (5 > 0, without loss of generality we set S = S = S and e = (5/2, and 
7>0,7<i^,7</i such that ( |10.16[ ) has a unique solution given by 



and the following estimates hold 



\xkit)\ < 



4k 



7v27r 



6 6^ 



-("'-DtRe-lfclp, 0<tR< 



(10.17) 



\Mt)\ < 



4k 



7V27r 



2C3C3 , 4^36-^ 



,(M-f)*H„-|fc|P 



-00 < < 0, 



(10.18) 



ll^llp-5.^-5 < 

dT 



8k 



(57V27r 
8k 



< 



C3C3 + -Cse 2 




C3C3 + ^C'se" 



Finally we have 



ll^llp-^,p-^-^ll(G"G)(0)t+ 



(57V27r 



C3C3 + ^Cge- 



1' 
SI ' 



ax 



p—6,p—6 



11/^ ^N/^MI 8k 

— ||(G-G)(0)||p + 



(527^2^ 



C3C3 + 76*36 2' 




Next we look at the solution of X ■ dqiY — p. Recall 



(3j = K (i?,(g') - B,) - [C{q') ■ d,,X{q') - C ■ dg,X{q')] - [c{q') -^-C-^ 



First we have 



a Tjl 



dsl 



dhl, 



.feGZ" -^J 
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and 



Substituting for the definition of Bj we obtain 



B. 



G -Ap— (5,cr- 



(10.19) 



Next, we write out an expression for C{q') ■ dq>X{q') — C ■ dqiX{q') given the following expressions 
X{q,t) = y{q) +T{q') G ApS.a-S, 

H\p', q, t) = G(p', q) + Tip', q') G ^p,,, 



We have completely solved for X{q,t). Recall from Theorem 9.1, Xk{t) is of exponential order with respect 
to tu- We differentiate and obtain 



and 



dq'X = (^dg'^X,...,dq'^^^x'j £ Ap-2S,a-25, 



/n+l 



n+l 



C{q') ■ dq,X = Ci.i{q')dq,X, Y Cn+lM)dq'X e Ap-2S,a-2S, 



\l=l 



1=1 



n+l 



{C{q')-dq,X)^ = Y.^,j{q')dq,X 



1=1 



n+l 



E 

n+l 

E 

i=i 

n+l 

E 

i=i 

5,;[:y(g) + r(g')] 



dq,[y{q)+T{q')] 



a2fV0 a2„l a2Ll 



E 

1=1 



7l!7(0)e*'' + -E7Sir(0)eU0e^'=- 



25 



+ 



+ 



7(0) + E o^^i^y'-'^^E ^^(0)4We-« 



n 

E 



feez" 



+ E a^r^(0)^^(*)^^'"^<..^('?') 



At this point we will pause to examine the terms obtained in the expression above. We begin with the first 
term 



E 

1=1 



f)2 frO 

— — (0) ^ ikme"'-'' 



E 



ii=i 



Jk-q 



The second term 



E 



- E ^(o)^"^^.'^(^) 



E 

1=1 



feez" 



i=i L fceZ" \mez 



E 

1=1 

E 



^ E ( E ^^(O)^'^^e- e-^ 



«E(E^^(o)(--H.y« 



The third term 

n 

E 



E S^(0)4(i)e^'="'9,;3^(9) 
kezn^Pi^Pj 



E 

(=1 



92 



E a^(o)^^(*)^^'" E ^^'y^^^'-' 



E 

1=1 L 

E 



E ( E ^Shio)eimyke 



i2?.l 



E (E^'^l7^(o)e-(*)(^-^)'y" 



.meZ" 



dp[dp^ 
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The fourth term 



E 



(=1 

The fifth term 



E 



■J feez" 



E 

feez" 



E 



E 



E 
E 

wez" 



feez" WeZ" 



52 „i 



E E^'^^^7^(o)("'-'^)'^— w 



The sixth term 

n 

E 



E 



= E 
= E 



E ( E i^^iOymit)kx,{t)e 
fcez" Vmez- ^^i^^J 



Jw-q 



The seventh term 
The eight term 



' dt 



,ik-q 



EE 



feez" Vmez 
The nineth term 



We finally write 



E'^'^^'^SrW + ^Ef E *7ST(0)(fc-m);yfe_, 



1=1 VmeZ" 
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and 



+ 



E 



" / o2 1 ' 

'^E E ^a^(o)(^-"^w-w 



+ «E E^a^(o)«-W(^-™W-'"W 



dv'n+idp'. dt 



-it) 



+ E 



dp'n+idPj 



1 dXk— 



AO 



ikWk p. , (0) 



+ 



= 1 \meZ" 



+ E 



^E( E z^(0)e^(i)(fc-H^y.-m 



-1 \mG2 



f pp. 1 ^ 

+ '^El E ^^(0)(fc-m),a.._„(i) 



i=l \m6Z" 



dp'n+idp'j dt 



meZ" 



dp'^^^dp'j dt 



-it) 



+ 



, dxk 



Ak-q 



Next we write out the expression {c{q') • ^ — C • . 

n+l 



"+1 "+1 /:)2f70 r)2/7^ 



=1 



= E 



dp'.dp' 



6 



and 



(10.20) 
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/'n+1 



/n+1 



, -' dp', dp'. 
We have thus separated the right hand side of 



(10.21) 



(10.22) 



into a quasiperiodic and a exponential-order-with-respect-to-time part. We first solve the quasiperiodic part. 
Assume 



Y,{q') = S,{q)+T,{q')^ S,{q) = ^ 5,,,e''=-^ T,{q') = ^ .F,,,W< 



fcGZ 



The j th component of ( 10.22 ) is given by 

\-d,,Y,{q')^Y.^^^'Ml')^Pj- 



n + l 



i=l 



From (10.19), (10. 2C) and (10.21) we obtain the quasiperiodic part of /3j 



n J: §^(0)e--+ E 
feeZ"\o -J' 



fceZ"\o Li=i 



Jk-q 



E (E«6|i%(0) 



Jk-q 



Zi\0 \;=i 



E 

feeZ"\o 



.|^(0) + E^%^M7(«) 



' dp'idpj 



n fl2 1 n+L „2 1 



E 



AO L 



m 7^ —k 

'^?7(0) + E^'^'2^''#4t(0) 



" 32 1 "+1 a2 1 



1=1 m£l 



Jh-q 



E 

fceZ"\o 



where we have set 



1=1 



1=1 me 



" dp'idp'^ 



n+1 



(10.23) 



1=1 



' dp'idp'j 



The quasiperiodic part of ( 10.22 ) is 



~X.d^{S,iq))^pf{q)= Yl P^^""' ^ A-2S. 

k£Z"\0 



(10.24) 
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Since we assume A G fir and since the right hand side of ( |10.24| ) has zero average, by lemma B.2, for some 
positive S < p , (10.24) admits a unique solution 



S,{q)= E ^fe^.e*" e ^p_35, 

feGZ" 



with iS," — and one has the estimates 



||5j||p-3d- < 7^ll/3f llp-2.5, 



dq 



< 



VJ 



p-3<5 



T5^ 



where tz7 = 2'*"+^^jJi±i)"^^. Next we want to solve the exponential-order- with-respect-to-time part of (10.22) 
From ( 10.19 ), ( 10.20| ) and ( 10.21 ) we obtain the exponential-time dependent part of Pj 



E 



fcG2 



-||(0)el(i) 



= 1 \m& 



1=1 



" / fl2 1 ' 



'^El E ^^^(0)eLW(fc-m)^x.._„.(i) 



= 1 Vmez 



kei 



Therefore, the exponential-order- with- respect-to-time dependent part of ( 10.22 ) is 



(10.25) 



To solve (10.25) we must find an exponential bound for the time dependent Fourier coefficients /3^(i). This 
exponential bound can be found by bounding each of the nine time dependent terms that make up Pfk{t). 
Note, to bound the terms involving sums over m we use the bounds found in appendix ^ For the first term 



|f(0)el(t) 



< ^Cie-('^-")*«e-l''^l'', < ci < ii? < oo, 





dh 



< ^C2e(^-^)*«e-l''^l'', -oo<tR<C2< 0. 





For the second term we refer to ( |10.15D 



^E( E ^Q^^{0)eUt)ik - m),y,^ 



1=1 \m£l 



30 



so that 



< 



WhlnWp I 1 



1=1 meZ" 



(52 



m\{p-S) 1-1 



meZ" 



meZ" 

^ E '^if (S)" l|G(0,g) -^(0)||,e-l'=-'"l(^-^^)||/i^||, |e^(t) 



„|fc— m|<5 



92 /ll 



= 1 VmeZ" 



^ E -;53f (S)"ll^?(0'9)-G(o)||pCi 

iGZ" 

'^^f(^) ||G(0,,)-G(0)||,C: 



g-|fc|(p-25)g-(^-|)tfi 



< ci < ifi < 00, 



«E( E 'd^.^^yrnmk-m),y,. 



1=1 VmeZ" 



1 K / n\n 



6(53 r 



< 



For the third term 



^ E 

meZ' 

3" K / n \ « / 1 



(^) ||G(0,g)-G(0)||,C2 



g-|m|pg-|fe-m|(^,-25)g-(^-|)tH 



r 

00 < tij < C2 < 0. 



-25 



||G(0,(z)-G(0)||pC2 



g-|fc|(p-25)g-(;^-|)tK 



< 



E ( eJ 



1=1 



o\k\S 



92jj-o 



so that 



" ^2^0 



1=1 



< 



E ( 



Jfc|<5 



dp'idp'j 



7(0) 



8k 



57 V 27r 



C3C3 + ^C3e-5-^ 




< tfl < 00, 



" 92JJO 

— (0) 
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< 



Jk\5 



8k 



J7V27r 







f.it^-i)tR^-\k\p 



-oo <tR<0. 



The fourth term 



meZ" 



so that 



«E( E *^(0)(fc-m);a=._W 



^ E 

meZ" 



8k 



< 



n „ / 1 1 1 1 p 



3"k 



GILA f 1 



C3C3 + ^Cae-'^^ 




g-\m.\p^-\k-m\{p-S) ^-{u-^)tR 



8k 



C3C3 + ^C3e-^^3 



e-|fe|(p-'5)e-(^-f)*K 



for < tR < 00, 



^E 



^ E 

meZ" 



/|1G||\ 8k 



< 



3„,ri|G|lp 



1 



G3C3 + 7^36- = ^^ 





-\m\p^-\k-m\ip-S)^{p.-^)tR 



8k 



G3C3+ 7^36-5^^ 




-mp-s)^-{fi-^)tR 



for —00 < ffl < 0. 
The fifth term 



'^E ( E ^^(0)e^(t)(fc - m),a.fe_„(t)^ 



< E - m||ei.(i)x,_„WI < E «^ (^) e^''-"'\'\\hlUeUt)xk-m{t)\ 



so that 



'^E(^E *^^(0)e^W(fc-m),Xfe_„(t)^ 
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g-|™|Pg-|fe-™l(p-'5)g-2(i'-|)tH 



< 



3" 



Ci 



8k 







g-|fc|(p-'5)g-2(i'-|)tH 



for < Ci < tfl < oo, 



^ E 



< 



C3C3 + -^C^e—--- 



g-|™|Pg-|fc-™l(p-5)g2(M-f)tH 



5^ V 1 - e 



-25 



8k 



C3C3 + ^Cae-s^^ 




g-|fc|(p-5)g-2(M-|)tR 



for —00 < tji < C2 < 0. 



For the sixth, seventh and eighth term we refer to lemma F.5 when deahng with the bound of the derivative 
of an exponential-time dependent function. The sixth term 



«2fV0 J 



< 



7(0) 



C3C3 + -Cae 



g-(i^-|)tng-|fc|p 



for < tii; < 00, 



OPn+lOp. dt 



< 



7(0) 



(5 Sj\/2tt 



Caca + ^Cae-^^^ 




g(M-|)tHg-|fclp 



for —00 < ifl < 0. 
The seventh term 



a^^^^ ^^^^ ^ 



< 



E- 



(52 



dt 



it) 



< 



E- 



IIGII 



(52 



dt 



-it) 



so that 



^ E 



< 



«l|G||p 
3"k||G||p 



(53 ,57^2^ 



CaCa + ^Cae-^^^ 




g-|»"|Pg-(''-|)tHg-|fc-m|p 



Caca + ^Cae-^^^' 
() 



g-|fe|(p-<5)g-(!^-f)tH 
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for < tfi < CO, 



dr>' dv' ^'^r^' 



< 



3"«||G||, 



C3C3 + 7^36-5"^ 





g-|'"|Pg(M-|)tBg-|fc-»n|p 



eS J 53 S^V2^ 



C3C3 + -,Cse—-^^ 




-|fe|(p-'5)g-(M-|)tK 



for —00 < tij < 0. 
The eighth term 



dxk 
Hi 



so that 



< 



meZ" 
rneZ" 



(57\/27r 



C3C3 + ^C3e- = ^^ 



g-(i/-|)tHg-lfc-mlp 



< 



3"kCi 



g(-|)<5/l\" 8k 



C3C3 + ^^36-^^^ 




Q-\m\P(,-\k-Ta\p^-2{u-^)tf 



5^ \eS J d-fV2^ 



C3C3 + jCse-^^' 




g-|fc|(p-«)g-2(.'-|)*fl 



for < Ci < fi{ < 00, 



^ dx 



meZ' 



< E 

meZ" 



g(M-|)«Bg-|™|p£ 



(5 SjV^n 



C3C3 + TC3e-^=^ 




= E 

meZ" 



kCi 



8k 



< 



3 t\jCrt 



5^ b-1\f^K 



C3C3 + 7^36- = ^^ 



g-|™|Pg-|'=-™|Pg2(/i-|)tK 



8k 



^3 \^e(5y 57\/27r 



C3C3 + ^C3e- = ^3 



g-|fcl(p-<S)g-2(,x-f)tB 



for —00 < fij < C2 < 0. 
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The nine th term 



"+1 ^)2f,l 



<«i?i 



(52 



so that 



^2 



Cle-("-5)*«e-l'=l^ < ci < iij < oo, 



< K|{j!l^fte''-'l'«e-l'l», -oo < (k < C2 < 0. 



Note we have bounded the sums on m as indicated in Appendix 3. Putting the estimates for the nine terms 
we obtain the following bound for (if . 



\Pfkit)\ < 







3" n f n \ 



||G(0, q) - G{0)\\pCi ] e-l'=l(''-2i)g-(<^-|)tR 



8k 







-|fc|(p-<5) 



3" 



1 



8k 



l-e-2V 



6(53 y I 1 _ g-25 



8k 



(57 V 27r 



(5 





g-|fc|(p--5)e-(i'-|)tB 



e-\k\{p-S)e-2ii^-^)tR 



7(0) 



e('^-i)'5 8k 
<5 (57-v/27r 



3"K|lGt ( ^ 



l\"e(^-i)* 8k 



3"kCi 



<53 6-fV2^ 
e(-l)^/l\" 8k 



Caca + ^Cae-^^^ 





GaCa + ^Cae-^^^ 



CaCa + ^Cae-5'=^ 





(5^ Ve^/ ^7\/27r 

g-Ci'-f )tfl^ < ci < tfl < oo 



-|fe|p 



e-|fcl(p-^)e-('^-t)*« 



e-|fcl(p-<>')e-2(''-t)*« 



+ Kl^llcie-I'^l" 

where we have used the fact e-2(''-l)*« < e-'^"-^^*^. We can further simplify by using the fact e-l'^l^ < 
e-mP-S) < e-|fe|(p-25) so that 



3" K / n \ » 



1 



6(5^ r\e(5/ \1 — e 
8k 



-25 



||G(0,(z)-G(0)||,Ci 



(57 V 27r 



Caca + 7Cae-5=^ 
(5 



1 



1-e- 



-25 



8k 



(57 V 27r 



GaCa + ^Cae-5^^ 
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similarly 



+ 



+ 



3" 



l-e 



-2S 



8k 



C3C3 + jCse-^^' 




t(0) 



C3C3 + ^C3e-5=^ 





C3C3 + lc3e-^'' 




+ 3"kCi- 



:(-|)'5/l\" 8k 



J3 



6(5 y SjV^ 







g-|fe|(p-25)g-(.-|)t«^ 0<Ci<iii<OO, 



8k 



(57\/27r 



C3C3 + jCse-^^' 




+ 



1 



l-e 



-25 



8k 



+ K 



+ 



3" 



l-e 



-25 



Co 



(57v27r 



C3C3 + ^C3e-5^^ 




'(57\/27r 



C3C3 + ^Cae-^'-' 




+ (^3"k||G||p 
+ \TkC2 



t(0) 



e(A'-|)'5 8k 



eS. 



5"^ (57V27f . 



C3C3 + -.C^e-i"^ 




C3C3 + ^C3e-5^^ 





8k 



5^ \eSj 5-fV2n 



C3C3 + 76-36-^^3 




g-|fc|(p-25)g(^-i)t„^ -00 < ifl < C2 < 0. 



At this point we must relate the constants used to bound the time dependence, Ci and C2, to a constant which 
bounds the entire perturbation. Recall that by assumption on the time dependence and by the quasiperiodic 
q dependence we obtain the following estimates 

\hl{0)el{t)\ < (7ie-l'=l^e-("-3)*«, < ci > < 00, 
\hl{0)eUt)\ < C2e-l'=l^e('^-i)*«, -00 < < C2 < 0, 



and clearly 



\hUO)el{t)\ < Cie-^'^-s)*", < ci > iij < 00, 
\hl{0)el{t)\ < C2e(^-5)*«, -00 < < C2 < 0. 
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Wc want this estimates also to inchide the interval C2 < Ir < ci. The assumption on the function e^{t) 
in this interval was simply to be a bounded function. Consequently, there exists a constant Eq such that 
Ci < Eq, C2 < Eq and the following estimates hold 

\hl{0)el{t)\<Eoe-\''\p, < < 00, 

\hl{0)ei{t)\<Eoe-^''^P, -00 < < 0. 

Now recall the total perturbation H^{0, q') — G{q) +T{q') and assume \G{q)\ < Eq- By the argument above 
we have 



6) 



feez" 



feez" 



Finally we have 



C\ Eq ^ Ej" ^ Eq -\- Ej" — E E\^ 



C2 ^ -^0 ^ Ej" ^ Eq -\- Ej" — E ^ E\ 



C3 =max(Ci,C2) < El, 
where Ex is the bound on the total perturbation H^{0,q'). The estimates on become 

1 " 



similarly 



\(3fk{t)\ < 



+ 



+ 



+ 



+ 



+ 



+ 



1 / 3" K / n 

^"^l^fv^y V 1-6-2.5 



1=1 

on I \lG\\p 



1 



1 - e 



-2(5 



||G(0,g)-G(0)||, 

C3 + ^e ^ 



C3 + 7e 




; V 1 - 



(57-\/27r 



ca + 76" 



7(0) 



l5 SjV^TT [ 



C3 + -e-5'==' 



1 \" ef'^-i)'' 



3"k||G||p — , ^„ ^ ^ 



2 



C3 + -e 



+ 1^1 



3"kCi 
1 



s(''-|)-5 / 1 \ " 



(52 



<53 Ve<5y ^7^277 
«£;^e-|fe|(^-25)e-(--i)*fi^ < ci < i^j < 00, 
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^( eS 



1=1 







3" 

|2 £V0 



1 



1 - 



(57-\/27r 



C3 + 7e 2' 





1 - 6-2-5 



Co 



5'f\/2TT 



C3 + 





92 if 



dp'n+idp'j 
3"k||G| 



7(0) 



" I 6(5 



5 (S7V27r 

' e(''-i)* 8 



C3 + tg"^"' 



3"kC2 



C3 + -6-3« 





63 + ^6 



"4 



(53 ^6(5/ (57%/2^ [ 

k£;i6^I'^I(''-2'^)6(^-^'*«, -OO < ifl < 62 < 0. 



We take a pause to a(idress the analyticity issue. We know from the functions that make up fifiq') that 
Pf{q') G Ap-2S,a-25- Our estimates above demonstrate the Fourier coefficients of (3f{q') depend on k as 
g"|fe|(p-2(5) .^yj^icJi agrees by Lemma A.l with our expectations. Next we want to further simphfy the estimates 
above by removing dependence on Ci, C2, G, k, H'^, and We first find an estimate for Recall from 
Lemma 5.4 there exists a positive constant / such that 



Also from ( 10.^ we have 



C 

A 



C 

A 



kEiCa 



>f 



>m\\ 



kB-C ■ dg'X 
kA 



< 4B\\p,a + \\C- dg,X\\p,, + k\\A\\p,, < 2nEi + m-^\\dy,X\\p, 



5 2'"=. + i3?£^l|G(0,,)-G(0)||, + -^ 



C3C3 + jC3e"'°" 



We finally obtain 



m < 



^^^^ + ^nrs^''^ 



8k 









C3H 


(5 





r7i(527v27r 

With this bound we can bound |Cl <(" + 1)11^11- Also note 

\\G{0,q)~^{0)\\p < 2\\G{Q,q)\\p<2\\H^\\p<2E,, \\G{0,q)\\p < E,. 
The bounds on ^, G together with kGi < kEi < 1, KC2 < kEi < 1, and 
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are used to reduce the bounds on to the foUowmg 

1 



Wk{t)\ < 



+ 



+ 



+ 



1 

3" 



3" 1 / 71 y 



6^T\edJ Vl 



1 



-25 



1 \ n 



eS J S S-fy^ 



C3 + Te 



1 

Vl- 6-2^7 ,^772^ 
1 



3" 

le^-^'" 8 r 2 _,,3' 





C3 + -e-. 
2 " 



C3 + tg" 



<5 (5 (57V27r 
3" ( 4 



e6 J S-fV2^ 
e(-l)Vl\" 8 



2 



C3 + -e 2=-' 



C3 + -:e 



1 + 
52 / 



2vj 











C3H 


(5 


(^1 





< 



3"16(n+ l)tn 
(54»+3r7/TO 



C3 + ^e-5« 




«:£;ie-l^l(''-25)e-(--|)*^ 



for < ci < tfl < 00. 

Now we use the fact for a; > 0, e^^ < a; + 1. Setting x = |c3 we have 

2 /<5 



C3 + ^e < C3 + ^ I ^S3 + l ) =2C3 + ^ < ^(C3 + 1) 



So that 



3"32(7i+ 1)ti7 



6 \2 



(C3 + 1) 



Siri+iYjfm ' (1 - e-*)^ 
for < ci < tfl < CX3. 

We further simphfy by using the fact S < 1 and the following 



S S 



I < I < 1 + 1 < 2 .-e 



The exponential bound of is 
for < Ci < t/j < 00 and 



(C3 + l)e'^Ati;ie-l'=l(^-2^)e-(''-i)*« = C''i«:£;ie-l'=l(''-2^)e-('^-i)*«, 



for — cx) < < C2 < 0. Next, by Theorem 9.1 there exists S > 0,6 > 0,6 > 0, without loss of generality we 
set 6 = 6 = 6 and e = (5/2 , and 7 < t^, 7 < /i such that ( 10.25 ) has a unique solution given by 
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Note, since 



J2 /3ffeWe*'e^p_2.>-2*, 



when using the bounds obtained in theorem 9.1 we replace S with 3S and the foUowing estimates hold 

2 



\^jkit)\ < 



8k 



bd 



5(57 V27r 



g-|fc|(p-25)g(M-|5)tH 



5S 



El, <tR< oo, 
El, -oo<tR<0, 



\\:f, 



j II p— 3(5,cr — 3(5 



< 



8k 



5(57v27r 
8k 

< 



55 



p— 35,fT — 3(5 



5527^2^ 



Si 



where C^^ = max(C'^'^ , C^^). Now that we have obtained estimates for Sj{q), Tj{q') and their derivatives, 
we can obtain estimates for Yj{q'), Y{q') and their derivatives. Recall Yj{q') = Sj{q) + J-j{q'). It follows 



\Yj\\p-3S,a-35 < J^II/3j^llp-25 



5S-/V2TT 



5o 



El 



and 



n7 Q 8k 

\\Y\\p^'i8,a-?.5 < ^^^IIPj llp-25 + 



dYi 



< 



p— 3(5,0" — 3(5 



5(57v 27r 



5o 



Si 



8k 



5(527V27r 



5o 



Si 



9y 



< 



ll/3fllp-2* + 



8k 



00 



Now we apply the estimates above to the generating function x = X{q') + ^ • g' + Yi{q')p[ . We have 



dp' 



n „Q|, 8K 



p— 3(5,cr— 3(5 



||>"||p-3(5,(T-3.5 < YpTiWl^J ll/='-2'5 + 



5(57v 27r 



C3 



5(5 



'3 „— o(5C3 



Si ( - ) . (10.26) 



Next we must relate the constant ||/3y^((z)||p-2(5 and Ei. To do this we find the following estimates for the 
terms that make up \\I3^ {q)\\p-2S- The first term 



/cGZ"\0 



^3 \feeZ"\o / ^1 



p'=0 



and 



a(G-G) 



< ^||G(0,(7)-G(0)||p-A- < ySi, yqeVp-25. 
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The second term 



feez^xo 1=1 



and 



" a2o-0 



fceZ"\o (=1 



< 



< 



2^ r9T){r9T)'. 



1 SpJSp; 



ri52n+l 



||G(0,g)-G(0)||, 



2wkEi 



y^ a^g" 



< 



2wK,Ei n 



where we have used in the last inequaUty Cauchy's inequaUty and the fact that ||g||p,CT < 1- The third term 
« E E f E ^IS7(0)*,»e-..-. - .||^(0).,») 



so 



^ i/ci^fce 



ik-q 



1^1 \kei 



" / ,q2„l /)2 1 



< 



UK 



\\y{q)\\,.s\\G{0,q)\\,.s + >^Yl W-kiP')\\p-S 



feez" 



<!!^^||G(0,g)-G(0)||,||G(0,<z)||,_, 



J2 r^2n I 



< 



|fe|(p-5) 



-|fe|(2p-5) 



The fourth term 



/"+! ^2„1 



E (E»6^(»)l»''' 



n+l 



n+l 

1=1 



d 



iGip',q)-G{p')) 



dp', dp 



p'=0 



{2p-S) 



fceZ"\o ' -J 



n+l 1 9 P 

<«El^'l^l|G(o,g)-G(o)||,<^|^|, 
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and with the four estimates above we obtain 

2k ^ 2wkEi n 2ntA7K^ 



2k2 \ f n\ 



Finally we have 



2kE 



9x 



dp' 



Next 



< 



p— 3(5,fT— 3(5 



2k 

T 



2zukEi n 2nzuK'^ 



2n-\^ny 



2kEi 



8k 



00 



1 , 



9x 



< 



p— 3(5,0"— 35 



ax 



We put together the bounds for ||^||, 
dx 



p— 35.17 — 35 
dX 



dq' 



+ {p-3S){n + l) 
and 



dY 



dq' 



2zuK ^ 



p— 35,(7—36 

8k 



dY 
dq' 



,.S5,a-35 r52n+l J2^^ 



2 ^< 



p— 35, a — 36 

4' 



C3 + 76 2 




+ ||C||p-3(5,o--3<5- 

p— 35, a — 35 

and obtain the following 

El 



8K(n+ 1) 



p^2n+l ll'^J 
2n7K 



5527^2^ 

8k 



55 



C3 + -e~^'' 




We set p = p ^ 3(5, ct = fj — 35 and examine the estimates for 
dx 



dp' 



and 



term of 



from (10.26), is less than the third term of 



dx 

dq' 



dx 
dq' 



4 

~5 

p,a 



E, 



Since S << 1, the first 



r52 



-\\Pfiq)\\p-25 <{n + l)-^\\pf{q)\\p_,s. 



Also, the second term of 
8k 



dp' 



00 



,from (10.2(:), is less than the fourth term of 
8K(n + 1) 



3p-i5c3 



El < 



00 



El. 



Consequently we obtain the following estimate 



dx 




dx 




9x 




dq' 


p,CT 


dp' 




dq' 


p,a 



2WK 



-E- 



8k 



/ n ^O/ Ml 8K(n + 1) 

+(- + i)t^I1/3^^('z)IIp-2. + ^^ 



r52n+l ,52^72^ 



5<5 



C3 + ^e- . 



^1 



2Ki;i 



2n7K 



8k 



r52"+i m527^/2^ 



C3 + ^e-^^^ 




El 



2VJK. „ 

< ^, Ei 



8k 



W2"+l ^ 



2 _i£3 

C3 + 2 



4\ " 



2k 2zukEi n 2nwK^ 2 

7- 1 + r52n+l ^ + rJ2«+2 1 
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+ 



^7 



2voK, 



(2p-5) 

8k 



8K(n+ 1) 
5(527^/27? . 



55 



O 77. OK 



C3 + 



-fC3 



2wK ^ 8k; 



2 _^£3 

C3 + ^e- . 



/4\" r2K 2wKEin 2ntuK;2 



+- 



+ 



8K(n+ 1) 



-^1 



2Ki;i 



(2p-5) 
2t<7K 



■1 H -El H 

'^mr52"+i '^m<527x/2^ 



C3 + TC 



1=3 



5527^/2^ 



00 



1 r„ ^ , 2zukEi 
2Kbj 



f 



El 
8kEi 



C3 + -^e 



|C3 



There are thirteen terms in 
the terms 



in the estimate above which we would like to reduce to one term. We first list all 

2WK 



1) 

2) 
3) 
4) 
5) 

6) 

7) 



rj2"+i 

8k 



^1, 



(527^2^ 

2tUK 
r^2n+2 

2w2k 



2 

C3 + ^e- ^ 



(n + l)i;i, 
■n(n + l)Ei, 



p2j4n+3 

2n72K2 
p2^4n+3 

r2eJ3n+3 



(2p-5) 



El 



9) 



10) 



4tzrK2 
52n+2r/ 

4^2^2 
mr2<54n+3j 

16tI7K2 



(n + l)X, 



(n + l)2i?2, 



m 



.(52"+4r/7V27f 



(n + 1)^ 



2 



C3 + ^e 



8K(n+ 1) 



5(527V^ 



C/33C3 + 1(7/336-1503 



55 



11) ^A, 

12) ^r/52n+l^l' 
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13) 



C3 + ^e 2' 



El. 



We begin to combine terms keeping in mind F, 7, 5, m, f, kEi < 1 and w » 1. The sum of terms 1 and 3 is 
bounded by 



2wK 



the sum of terms 2 and 13 is bounded by 

16k 



C3 + -,6-^"' 




the sum terms 1, 3, and 4, knowing that vcj > 1 and assuming without loss of generahty that F < 1, is 
bounded by 



the sum of terms 11 and 12 is bounded by 



{71"^ + 2n + 2)Ei, 



4WK 

mr/(52»+i ^' 
the sum of terms 1, 3, 4, 11, and 12 is bounded by 

the sum of terms 7 and 8 is bounded by 



{n+iyEt, 



the sum of terms 7, 8, and 5 is bounded by 

2 2 

-^,^^(l0n^ + 18n + 8)Ef 

the sum of terms 7, 8, 5, and 6 is bounded by 

1 



1 



1 _ e-(2p-'5) J m/r254"+3 



(10n^ + 18n + 8)S/, 



where, to bound 6, we have used the fact | (^) < zu. The sum of terms 1, 3, 4, 11, 12, 7, 8, 5, and 6 is 
bounded by 



1 



1 - e-(2p-'5) J TO/r2J4n+3 

the sum of terms 9, 2, and 13 is bounded by 



(lOn^ + 18n + 8)fi;£;i, 



(10.27) 



16ro(n2 + 2n + 2) 



/mr7,52"+4V2^ 



and the sum of terms 9, 2, 13, and 10 is bounded by 



^ -J tul6(n^ + 2n + 2) 



/mr7(52»+4v^ 



C3 + tc" 



(10.28) 
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We use the fact 



2 2, 



together with 



1 



1 2 

1 - e-(^p-^) ^ 2p-6'^-'p-S'^-J^^^^]^ 



and combining ( |10.27D and ( |10.28D we obtain 



16w2(10n2 + 18n + 8)2 , ,^/4^ ^ 



Recall 



Then 



Finally 



where 



3"2i2n73(6n + 6)3(c3 + 1)^ 



6n + 6 



A 



H = lOn + 9, 

_ 3"22"+i2n73(6n + 6)^(c3 + l)' 



A 



(10.29) 
(10.30) 
(10.31) 



With the estimates on the derivatives of the generating function one can now estimate the analyticity domain 
of the flow associated with the corresponding Hamiltonian vector field. We have obtain the following 



6n + 6 



where 



p = p — 3S, a = a — 3S. 



In particular, note by (10.1) it follows rj < 1 since it can be written in the form 

~ nil. 



and therefore 



4(n + l)' 



6n + 6 



Consequently, Lemma C.2 applies and x generates a canonical transformation 



(10.32) 
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where p' = p—5 = p— 45 , and a' = a — S = a—4S. We have seen X G Apss.a-sSi Yi G ■Ap-j,5^a-?,s, for i = 
1, n + 1 and since X = <l' + Yi{q')p'i it follows x G •Ap-z5,ij-35- Next, we consider the following 

estimates. 



II^'IIp'..' < 1 



From the second inequality of Lemma 0.2 it follows |l-ff'|lp',o-' < ||i?|lp.o- < 1 



Fe^lp,,, \\UF-F\\p,,^, <r^\\F\\ 



From Lemma |C.2| we have the following estimates 



\\UF - F\\p^s,a^s < 4(n + l)^||F||p,,, < \\F\\p^,, 



Using ( lO.Sj ) we have 



and combining these estimates we obtain 

\\UF-F\\p,,,, <7m\p,^<v\\F\\p, 



(10.33) 



m! = m — 



The approach is to construct m' so the following inequality hold ||C'w||p,o- < m' Vu G C". By 

definition it follows 



Applying Cauchy's inequality, ( 10.33D and the fact that a' > we obtain 



r'2 



from which it follows 



Using C" = C + (C - C) we have 

||C'«|lp,,., < \\Cv\\p,,,, + UC'-C)v\\p,,,,. 
Using the estimate ||Cti||p',(j' < ||Cti||p_cr < "'^ll'f^ll, along with ( |l0.34 ) and ( 10.35| ) gives 



\\C'v\\p,.<im-^ + '-^^)M. 

Applying the inequality + b < {a — b)^^ for < b < a < 1, the estimate becomes 

4(n+l)?A"\ 



\\C'v\\p,^,, < 
This estimate provides the value for m' 



< m — 



(10.34) 



(10.35) 



4(n + 1)?7 
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/' 



We begin with the known condition 



Furthermore we have by definition 



C X' 
A 



> /IK'll 



(10.36) 



ciM) ~ aAi') = Q^i^H H]{o,q'). 



Applying Cauchy's inequality, using (10.33) and the fact cr' > cr*, gives 

\\UH-H\\p,^,, ?7 



C"i,7 ~ Ci 



< 



at ai 



from which it follows 



(10.37) 



Also one has A' = {1 + kEi(^a)X. From this follows A' — A = kEi(aX, and we obtain the estimate ||(A' — A)u|| < 
We now find an estimate for (a- We first find a bound for 

kEiCa 



We have 



and 



kEi(a 



c \^ \ V KB-C-d„,X 



A 



nA 



kB-C ■ dq'X 

kA 



< k\\B\\p,, + \\C ■ dg'X\\p.25,.-2S + n\\A\\p,^ 

< 2KEi+m-^\\dg,X\\p^2S.a-25 



^ 2«£;r + -^||(G-G)(0)||,+ 



< 2kEi 



■ujkEi 



8k 



C3 + -^e- ^ 



2 _i£3 

C3C3 + -Cae 2 




We know 

Its not hard to show 
Therefore 

^EiWCaW < 



kEiCa 



< 



C A^ 
A 

d A^ 
A 



>f\\v\\ 



+ 



2 



C3 + ^e- 



47 



< 



4"lln7 



Using ( |10.37D and ( |10.38D we obtain the following 
A- A' 



(C3 + 1)kEi = TkEi. 



V2 



{C - C")vi + (A^ - X'^)v2 
(A - X')v2 



(10.38) 



< 



[C - C) wi II + II (A^ - A'^) W2 II + II (A - A') ^2 II 



< (!l±i)2|K|| + i 



8 


2 5c3 ■ 






C3 + -e-^ 


m 



kE'i 11^21 



2 _Sc 



C3 + ^e 



k£;i||u||, (10.39) 



where u = (wi,W2) and clearly ||wi|| < ||u|| and ||w2|| < H^'H- Using ( 10.39| ) and ( 10.3(]| ) we obtain 



C A^ 
A 



A - A' 



> 



/- 



[n + 1)?7 



/ /mr(52n+i 



2tu \ IGk-Ei 



m5'^^\/2TT 





2 ^03' 






C3 H 




m 


\\v 



\\v\\. (10.40) 



For any matrix M one can write M' = M + (Af - M) and ||M'|| > || Af|| - ||M' - A/|| . Together with ( |10.40D 
we can bound from below as follows 



a A'^ 

A' 



> 



{n + 1)77 



2vj 



16 



> {f-^-^-v-v]M>[f- 



mS^j\/2Tr 
{n + 3)77 



C3 + 2 



lkll = /'IK'll 



With the assumption max(|| A||p,cr, ||-B||p,cr) < Ei, the following estimate holds 

||P||p^<, ^k\\A + B- p'Wp.^ < (n + 2)kEi. 
We want for the new perturbation P' a similar expression namely 

llP'llp,,,, =K'P' + B'-p'||p,,,, < {n + 2)K'E[. 



kEi 



(10.41) 
(10.42) 



max(P'||p.,,.,||B'||p.,,0<£^l 



To arrive at this estimate we must first find estimates for k'||A'||p'^o-', K'||i?'||p'^CT' and identify the values 
for k' and E[. We know 

K'A'iq') ^ H'{0, q') -a' = H'{0, q') - F(0), 
H\p', q') = U + P + {x, U} + [{x, P}+UH-H- {x, H}]. 
By our choice of x we obtained 



H'{p', q')^a+{l + eCA)X ■ P' + O(llp'lr) + [{X, P}+UH-H- {x, H}], 
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n'A'iq') = H'{0,q')-a' 



= a + UH-H-{x, H}] (0, q')-a- [{x, P}+UH-H- {x, if }](0) 

= ({x, P}+UH-H- {x, H}) (0, q') - ({x,P}+UH-H-{x,H}) (0). 

Therefore wc obtain the estimate «;'m'||p/,o-' < 2||{x,-P} + UH — H — {x,H}\\p/^„i. Using the previously 
obtained estimates 



\\{x,f}\\p-s,a-s<2{n+l) 



\p,crt 



\\Uf - f - {x,mp-5,a-6 <2.2{n+lf 



^ I \\J \\P,ci 

\\P\\p,a < kII^IIp,* + (n + <{n + 2)kEu 

p' = p-5 = p-A5, 



45, 



5 -6n + 6' 
we obtain the following 

IKX, P}+UH-H- {x, H}\1,,,, < IKx, P}||p',.' + \\UH -H-{x, H}]],,,^, 
< \\{X, P}\\p-5,a-S + \\UH -H-{x, H}\\p_s,s_s 



< 2(n + 1) ( ^ ) ||P||^,, + 32(n + 1)^ f ^ ) \\H\\~,,, 



< 2(n + 1) ( ^ ) (n + 2)kEx + 32(n + 1)^ ( ^^^^ 



< 2(n + 2)\Er ( ^ ) + 32(n +lfi^ 



< [l + 32(n+l)2] 



6n + 6 



< 



l + 32(n+l) 
(6n + 6)2 



2 1 



2 ^ 2 



Where we have used the 



We then have 



2(n + 2fKEx < 



6n + 6 



(10.43) 



(10.44) 
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Next recall 



Using Cauchy's inequality we obtain 



^'\\B'\\p'^.' < ^m' - H'mwp',^' = ^m' - amw^,,,, - 1\\{^^p}+uh -h- {x,h}\\,,,,, 

a a a 



? 2 



A 



(10.45) 



Finally with (|l0.44D and (|l0.45|) we set 



and clearly max(|| A'||p/,<,-, \\B'\\p,^„,) < E[. 



Every time the iterative lemma is applied an equation of the fo rm o f ( 10.8 ) must be solved. The C,a that 
is chosen at the first step must satisfy the equation of the form of ( 10.8 ) for all other steps. Therefore at the 
step fc + 1 we have 



A'^+i = (1 + n''EicA)\'' = (1 + «'=i?i'a)(i + a-)A° - (1 + a+i)A^ 



(10.46) 



Therefore we need an estimate for 1 + Cfc+i- 



i + C 



Assuming i < 1 + ( < |, with (10.38) and 1 + ^' (1 + kEiCa){1 + C) it follows 



1 + C - i^EiQAil + C)<1 + C'<1 + C + «;^^iCa(1 + 0- 
As will be shown later,we can assume k to be small enough so that 

i < 1 + C - \^^E,T < 1 + C' < 1 + C + \^E^^ < I 



where T is as in ( 10.38| ) 



L' 



We know A' = (1 + C')A°- Then it follows 



|A'| = |(l + C')A"|<xio = i'. 



11 Conclusion of the Proof 

We now describe the set up for several steps in the application of the iterative lemma. 



Zeroth Step 
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We begin with the given positive constants 



r, po,o-o,mo,«;o < 1, P* < Po, cr* < cfq m^^Kmo, Lq, E° 



and 



-\2mo - Lo\ — 

A (Jn 



We assume a Hamiltonian of the form described above defined on D 



Ho{p',q') = Uo{p',q') + Po{p',q') 

Uoip', q')=a^+X'-p' + lYl + R'iP'^ l') 

Poip', q') = i^oA%') + KoJ2 BtWi 
where all the functions are in ^pq.ctoi € fir, with |Ao| < Lq and the following bounds hold 



A ^ 

ma^{\\AX„,„\\BX„,,) < E°. 



We apply the iterative lemma. 



First Step 



Applying the iterative lemma for any (5o > such that 

Po — ASq > p*, (To — 4(5o > a* , 

and kqEI small enough that 

mo - 4(n + l)rjo/al > m., /o - ^^Hl^ > 0, 
there exists a canonical transformation 

(f'l '■ -0^0-4(50, CTo-4(5o ~^ -^po,<^o' 'Pl G ^po-4i5o,cto-4(5o 

that transforms the Hamiltonian into 

= iJo o 01 = Uiip',q') + Piip'.q'), 
which can be decomposed as done previously 

C/i(p', q')=a' + X'-p' + \Yl ^Ul')PiP'j + R\P', q') 
P^{p' ,q') = K^A\q') + K^Y,Bl{q')p\ 
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with the same T and with the new constants 



Pi = Po - 4(5o > p*, (Ti = (To - 4(5o > (T*, mi = toq - 4(n + 1)^ > m,, /i = /o 



{n + 3)?7o 



r 3^ A 



A 



kqE^, Ki — — TsW^cn ^1 



with ah functions analytic on Up^.o-i, and bounds 



A (Elf 



fiM < 



A 



Vw G C 



n+l 



<mr'|k||, V«eC"+i 
max(||Ai||p,,,,,||Bi||p,,^J < El 



From the iterative lemma we know the following estimate holds 

1 - ToAto^^l" < 1 + Ci = 1 + koE°Cao < 1 + ToKoE^. 



Later we will show that in fact 



i < 1 - ToKoE° < 1 + Ci < 1 + ToKoE^, < ^. 



Using ( |10.4et ) and ( |ll.3D we have 



|Ai| = 1(1 + Ci)Ao| < (1 + ICiDIAol < = i,. 



Moreover, for / e Ap^^^t, we have \\Uif - fWp^.ai < ?7o||/||po,cro, where Uif ^ f o 



Second Step 



Applying the iterative lemma for any (5i > such that 

pi - 4(5i > p*, (Ti - 4(5i > CT*, 

and KiE\ so small that 

mi - 4(n + 1)77i/ctJ > m,, /i 



(n + 3)771 



>0, 



there exists a canonical transformation 

02 : Dp, 

— 4i5i ,(Ti — 4(5i — 4(5i ,(Ti— 4(5i 

that transforms the Hamiltonian into 

H2=HiOcb2 = U2{p\ q') + P2{p', q'), 



(11.1) 
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which can be decomposed as done previously 



P2{p\ q') = K2A^{q') + B-{q')p- 
with the same F and with the new constants 

fj _|_ 3^7y 

P2 = Pi - 4(5i > p,, (72 = cTi - 4(5i > cr*, m2 = TOi - 4(n + 1)— > m,, /2 = /i ^ , 

with all functions analytic on Dp^ ^a-i , and bounds 



/2||f|| < 



C2 



A J 



max(|lA2|lp,,,,,|li32|lp,,,J<i?2. 



From the iterative lemma we know the following estimate holds 

1 - koE°Cao ~ ^CimEl < 1 + C2 < 1 + hoE°iCao + JCi'«i£^i • 



Later we will show that in fact 



i < 1 - «oS?Ca - Ir^KiEl < 1 + C2 < 1 + hoE^iCa + ^TiKi^i < ^. 



(11.2) 



Using ( |10.4eD and ( pT3| ) we have 



IA2I = |(l + C2)Ao| < (l + |C2|)|Ao| < -io-i* 



Moreover, for / G Ap2,a2, we have ||ZY2/ - f\\p2.a2 < m I1/I1pi,(ti , where ^^2/ = / o 02- 



(fc + 1)*^ Step. 



We repeatedly apply the iterative lemma k + 1 times. We choose 5k at each step small enough so that 



Pk - 4:Sk > p*, <Jk- 4(5/c > cr*, 



and KkEi small enough that 



mk - 4(n + l)rik/al > to*. 



jn + 3)?7fc 
/fc -5 > 0. 



53 



There exists a canonical transformation 

0fe+l : Dpk-4:Sk,cyk-iSk Dp^^cTk, 4>k+l e Ap^-4:Sk,<Tk-4:5k, 

that transforms the Hamiltonian into 

Hk+i ^ Hko (pk+i = C/fc+i(p',g') + Pfe+i(p',g')) 

which can be decomposed as done previously 



with the same V and with the new constants 



Pk+i = Pk - 44 > p*, (Tfe+i = o-fc - 44 > o-*, mfc+i = mfe - 4(n+ 1)— > m* 



/fc+l = fk 



(n + 3)7/fe 



ifc+i — t;^o, Vk+i — — ' 



A 



Kk+l 



k+1 



A (£;[■) 



fc^2 



with all fmictions analytic on D p^,^^^„^^^, 



' f>l6l a* 
with bounds 



fk+i\\v\\ < 



C<k+i 
X 



max 



Pfc + l,o"fc+l ' 



k+1 



From the iterative lemma we know the following estimate holds 

3 3 3 3 

I-kqE^Cao ^Ck-iKk-iE'['^ ~ -CkKkE^ < 1 + Cfe+i < l + KoE^CA„ + --- + ^Ck-iKk-iE'l-^ + -CkKkE^. 

Later we will show that in fact 

i < 1 - koE"^Cao ^Tk^ikE'^ < 1 + Cfc+i < 1 + ^^oE^Cao + ■■■ + ^TkKkE^ < | (11.3) 



Using ( |l0.4e| ) and ( |ll.3D we have 



|Afc+i| = |(l + a+i)Ao| < (l + |a+i|)|Ao| < -i* 



Moreover, for / e Ap^^,,a^^,, we have \\Uk+if - /Wp^+^.a^^, < ^ll/llp^.a,, where Uk+if = f o cj)k+i. 



54 



Limit as fc ^ oo 



In the limit case we formally obtain the Hamiltonian 
which can be decomposed as before 

Poo{p', q') = eooA°°{q') + e^Yl ^TWi, 

with constants A°°, p^o, ctoo, moo, foo, -^00, ^^oo) E'f" where all functions are analytic on D^^^,^^ 
with bounds 



/oolkll < 



A 



max(||A~||,„,.^,||S-||,„,.„)<i;r>- 

For the expressions in the limit fc — > 00 to have meaning it is necessary to show 

3 

Poo ^ P*i Coo > Cr*, TTloo > 'Tl*) Loo — l^^^ ~ ^ ^* 

which will allow to specify p*, cr* and m*, and will lead to Koo = 0, i.e. the perturbation vanishes in the 
limit. 

It is also necessary to show that the sequence of canonical transformations defined by 
where 

(t>k ■ Dp^,a-k Dpk-i,crk-ii 

converges to an analytic, canonical transformation. 

We have seen the following recursion formulas amongst the different constants 

AX AX M I f f (« + 3)77/0 
Pk+i = Pk- 4()fc, (Tfc+i = cTfc - 4:dk, mfe+i =mk- 4(n + Ij-o , Jfe+i = Jfe 5 

Lk^,-^L„ ^k+,-j^^n„ E, nk^.E, --. (11.4) 

We also have Kk, and r]k related by 

% = (11.5) 

We can use the expression for Kk+iE^~^^ to eliminate % and thus obtain the following 

Pk+i ^ Pk~ 44, (11-6) 
o-fc+i = (Tfc - 4^fe, (11.7) 
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TOfc+1 = ruk 



fk+l — fk 



3/2 



(n + 3) 

3/2 



fc+l\l/2 



— 2-^0, 



fe+i _ 1 J^kHkEi 



Relations ( 11. 6| ) through (11.10) can be summed to give 



Poo = Po - 4 4 > P* , 



fe=o 



O^oo =0-0 — 4 ^ > cr* , 



fe=0 



_ (n + 3) ^ / .. /. 

Joo — Jo 5 2^ ('^fe+l-C'i j > /* 



fc=0 



moo = "^o 



4(n+ 1) 



k=0 



Note that using ( 11.4 ) we can rewrite ( 11.5 ) as follows 



r2f< 
Oh 



(11.8) 
(11.9) 
(11.10) 
(11.11) 

(11.12) 

(11.13) 
(11.14) 
(11.15) 

(11.16) 
(11.17) 



Defining the sequence {nkEi} serves to define the sequence {6k}- Now we must make an appropriate choice 
of the sequence {sk}- 

First we assume the sequence {KfeiJf } has the form 

KkE'l^CkKoEl A: = 0,1, 2... (Co = 1). 

The choice of sequence {wfeijf } is important since it will result in bounds for EqEi which is the size of the 
initial perturbation. We define the following series 



f cl 



2 \ 2 



k=0 k=0 
oc oo 



\Ck+i 



(11.18) 

(11.19) 



k=0 k=0 

With the above notation ( [11.12| ) can be rewritten as 



po - P* > 4 



A2 



(«;oi??)'/'^s or KoE°< 



2N 4 „4 

Po - P* \ mli^a^ 



where we have used the fact m* < ruk (which gives a more stringent bound). 



A2 



(11.20) 
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With the above notation ( [11.13| ) can be rewritten as 



(To ^ o"* > 4 



A2 



(«;oi?i°)'/'^s or KoE°< 



\ 2ii 4 j-4 



4s 



A2 



Similarly (11.16) can be rewritten as 



or 



and (|11.15| ) as 



(Jq (mp - m^)^ 
16 (n + 1)2^2 ' 



/o - /* > 



^^(.0i??)V2t> (!i±^(.0i??)l/2i 



^0^ 



(n + 3)2 



/o ~ /* 
t 



For the sequence {C^} we make the choice 



so (11.18) and (11.19) become 



fc=0 



= 2^2-2fe = 



fc=0 



oo ^ 



fc=0 



1 



< 1. 



2-N y 2^-1 

We therefore obtain fom- estimates for eo from ( |11.2C| ) ( |ll.2l| ), ( |ll.22[ ), and ( |11.25D 



As J A2 

2N 4 „4 



E^O ^ / CTo - cr, \ mj^a^ 



V 4s 
crj^ (mp - m^)2 
''''^1 ^ l6 (n+ 1)2^2 ^ 



A2 



(n + 3)2 V < 



/o - /* 



We now make the following choices 

Po 

P* = TT: 



TO* = — , /* = y and cr* = 



so (11.28) and (11.27) become 



»o^f<(|) 



32/ 28 A2 ' 
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Since ctq < po 



Also we see 



.327 28 A2 {32) 28 A2 



(To (too - TO*)2 (Tq TOq 



16 (n + l)2t2 26(n+l)2' + 3)2 \^ t 



/o — /* 



4(n + 3) 



f2 

Jo I 



and is not hard to see 



32/ 28 A2 26(n + l)2' 



2*^ TO;^/oV, ^0 r2 

V32/ 28 A2 4(n + 3)2''"' 



Therefore ( 11.31 ) gives the maximum upper bound for the size of the perturbation, kqE^. 



Now we develop the estimates for C,k+i with the assumption at the zeroth step —1/2 < (^a < 1/2 and 
then take the limit as fc — > cx). 

Assuming at the zeroth step —1/2 < (^Aq — (a < 1/2, at the first step we obtain 



1 + Ci < 1 + koE°Cao 



At the second step we obtain 



I + C2 < I + Ci + ^CaiKi^^i < 1 + koE°Cao + ^Ca^kiEI. 



Continuing in this manner we obtain in the (fc + 1) step 



1 + a+i < i + Ko£^?ao + ^^CA.«;«-Bi. 



i=l 

One obtains a similar expression for the lower bound and combining it with the upper bound we have the 
following 

1 - hqE'ICao -^K<1 + Ck+1 < 1 + ^oE"^Uo + l'^, (11-33) 

where k = X^iLi CAiHiE\. 

Now we look for a condition on k. Using the estimate on \C,a\ we have 

l + noE'l\Uo\ + \\n < l + «oi?i"To + ^|«| <l + ^£. 



where e — X]i=o '^i'^i^i- Similarly for the lower bound 



l-«oi??lU|-^N > l-«oi??To-^|«|>l-^£. 



Therefore expression (11.33) becomes 



so that 



We want the following to hold 



3 3 
1 - -e < 1 + C/c+i < 1 + 2 



ICfe+il < 



3 1 
2"<2- 



(11.34) 
(11.35) 
(11.36) 
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Therefore we want e to be a converging series with the above upper bound. 
With the previous definitions of KfcE'j"', Ck, the foUowing estimates hold 

hm -e ^ -} Tk^kEl < - > — rv-^ 



k=0 



Using ( |11.17| ) in the estimate above we obtain 

3 / 4"llti7 



ICo 



< 



< 



2 \r7m*/. 



m 



oo 

E 

k=Q 
2 f2 



3 / 4"lln7 ^^ mlf§y/^ 
2 vr7TO*/* 



A 



OO 



/2 



/£ = 



3 / A"Uu7 \ mlf§^ 
2 \r7m*/* 



A 



We therefore require 



3 / 4"lln7 ^^ mlfgy/a 
2 [r^mj. 



A 



0\l/2 



< 



2' 



and we obtain another restriction on the perturbation 



1 /rw* 

° ^ 9 V 4"lltJ7 



A 



(11.37) 



We now compare (11.37) and (11.31) and choose the smallest of these bounds to be the bound on kqE^. 
After substituting for A in (11.37) we obtain 



2^^w^6n + 6)«(C3 + l)4j;2('-+M) 



and after substituting A, ( 10.31 ), in ( 11.31 ) we obtain 



32 



28 32«24"+24n76(6n + 6)^(03 + I)4e2('^+M) ' 



(11.38) 



(11.39) 



A simple inspection reveals (11.3E) is smaller and therefore the upper bound on the perturbation kqEi. With 
( 11.3!^ we obtain the bound 



ICo 



< 



< 



3 ^^4-llzn ^ ^§/^V^(^^^0)i/2 



2 

32 



r 

/ 3"2i5n72 



(6n + 6)4(c3 + l)2e''+*'A' 



12 Time Aperiodic Perturbation Tending to a Quasi-periodic 
Time Perturbation 

In this section we consider a nearly integrable system generated by a Hamiltonian consisting of an integrable 
part and a small perturbation. As before we consider a real valued nearly-integrable Hamiltonian in action- 
angle variables of the form 

Hip,q,t)^H''ip) + KH'ip,q,t), 
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where p = (pi, ...,Pn) G -B C M", q = {qi, ■■■,qn) G T"^ arc. respectively, the action and angle variables, and 
K € M is a small perturbation parameter. The form of the perturbation is as follows 

H\p,q,t) = Q{p,q,t) + E{p,q,t), 

Q{p,q,t)= 5fee*(«'''*\ 

E{p,q,t)= J2 /fe(p)efc(i)e*«, 
feez" 

where Q{p,q,t) is time quasi-periodic , with = {6i,...,6m) as the vector of basic frequencies. f{p) is a 
bounded function and ek{t) decays exponentially to zero as time goes to infinity. Consequently, as time goes 
to infinity the entire perturbation tends to a t-quasi-periodic function. 

The result of the theorem will be similar to the KAM-type result obtained previously in this paper. 
Mainly, we will prove the preservation of cylinders of the form T" x R where the tori T" can be identified 
with the tori of the integrable system generated by H^{p). We transform the non-autonomous system into an 
autonomous system making each of 6it, 6mt a dependent variable qn+i, qn+m respectively. We will use 
either notation throughout depending on need to clarify. We add conjugate variables r = (n, r™) which 
are needed to obtain a Hamiltonian form. We will sometimes use the notation p' = {p,t) and q' = {q,9t). 
The Hamiltonian takes the form 

H{p,q,t) = e-T + H\p) + nH\p,q,t)=H\p,T) + KH\p,q,t), 

and Hamilton's equations are given by 

^ _ dH _ dH^ 



dq dq 
dH dH^ 



d{et) d{9t) ' 

^ dp dp dp ' 



The angular frequencies are given by A(p') = ^^ip') = {^^, ■-, ^-,0i, ...,6m)- We will often use the 

dpi 



notation A = (A, 6*) where A = (|^, |^). We define the complex extension of M" x x T" x M"' as 



follows 

Dp,a,po = {(P, T,P, et) G C^^+^^l \\p - Poll < p, \\t\\ < a, Re g e M" mod 27r, ||Im q\\ < p, |Im t\ < a}, 

and define Ap.a.po as the set of all complex continuous functions functions defined on Dp^„,p^, analytic in the 
interior of Dp^^^^p^ and real for real values of the variables. The normal form is again given as the Taylor 
expansion of the Hamiltonian about p' = 0. Namely 



_j n+m 

H{p, T,q,et) = a + X-p' + A{q') + B{q') -p' + ^Yl ^^'O WiP'j + W , q'), (12-1) 



where 



a = H{0)=HO{0) + KH\0)=H%0) + Km{0), 
KA{q') = H{0,q) - a = K{H\O,q,0t) -B^{0)), 
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and 



aeR,A, B„ a,,,R e A^,,, R = 0(|b'11). 
The matrix C* will be defined as before 

and the same non-degeneracy condition will hold 



det 



dpjdp. 



(0) ^ 0, 



ord||z;|| < ||C'*w|| <d ^\\v\\, Vw G C" , for some positive constant d. As before, this non-degeneracy condition 
will imply an isoenergetic non-degeneracy condition as proven in lemma |3.3[. In other words, although the 




matrix C* j — gpi^pi (0) i, j = 1, n, n -I- 1, n + to is a bigger matrix, the matrix I defined as 



C* 

X 



can be shown to be nonsingular and therefore a lower bound can be found for ||Iti||p.cr, \/v E ([^n+rn+i^ Using 
this lower bound, a similar constant, /, can be found to bound the expression 

A 



as was done in lemma 5.4. Ultimately we will need the implied invertibility of the matrix above to solve 
completely for the generating function. 

Because the form of the Hamiltonian perturbation in this section is different from the one consider 
previously in the paper, we will re-define what we mean by p', g'-exponential form. In this section, functions 
of p' , g'-exponential form will be understood to have the following form 



Fip',q') = fip')+rip',q')+g{p',q'), 



r{p',q')= Skip')e''-^'''''\ 
9{p',q')= /^febOefcWe*', 

where F{p',q') E Ap^a and as before ek{t) is of exponential order with respect to time. Clearly the only 
difference in the definition is the quasi- periodic time dependence added to the term r{p' , q'). Before examining 
how this modification will affect the lemmas proven in section ^ we prove a lemma concerning the new form 
ofr{p\q'). 

Lemma 12.1 

Given a function r(p' , q') E Ap,a of the form 

rip',q')= Y Skip')e^'<^-''\ 



61 



it follows r = so(p')- 

Proof: We define the vector k = (/ii, /12) G Z"+'" where hi G Z" a^ii /12 G Z™ and 



■1-1 -I L + " 

1 



dqdt 



dt 



= so(y) + 



E 

/i2ez™\o 



1 /-^ 



„ih2-9t 



dt 



S(o,,,)(p') ^lim^ ^ = so(p') □ 



/i2eZ"\o 



2T 



The result of lemma [F.l| still holds. With the result of lemma |12.1| , lemma F.I still holds. That is, given 



it follows 



^ /ife(0)efc(t)e^^« 



The results and proofs of lemmas comparable to lemmas |F.3| and FA do not change and follow the same 
arguments. A lemma comparable to lemma |F.5| concerning the time derivative of a function F(j! , q') of 
p' , (/'-exponential form would follow the same line of arguments but one should take in consideration the 
time derivative of the function r{p',q') = Sfe(p')e*'^'^'''^*^. Since such time derivative results in a g 

periodic time quasi-periodic function, the total form of dF{p' , q')/dt will be of p', g'-exponential form and the 
domain of analyticity would be as stated in lemma |F.5| . Finally, the added time quasi-periodic dependence of 
r(p', q') will not change the result of lemma F.7 concerning the p' , g'-exponential form for the Poisson bracket 
{x, F}. The arguments of the proof would follow the same and some sums would have to change from X]j=i 



to account for time derivatives that in lemma F.7 are zero but now must have to be added. 

The statement and proof of a lemma comparable to the iterative lemma [10. 11 would follow exactly the 



same up to the generating function used and the definition of the perturbation (10.12). The generating 
function will have almost the same form 

X^Xiq') + i-q' + Yiq')-p', 



X{q')=y{q')+T{q'), 



{q,6t) 



Y,iq') = S,iq')+T,iq'), 

S,{q')= J2 :F,{q')=Y.^>^^^(iy''' 



where now y{q') and Sj{q') depend quasi-periodically on time. The new perturbation must have the following 
form 



H\p',q')^G{p',q') + T{p',q')eAp,a, 
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where e\{t) will be assumed to have the same properties as before. Solving for X{q') involves solving the 
equation X-dq'X(q') = k (^A(q') — , which can be separated into a quasi-periodic part and an exponential- 
order- with- respect-to-time part. These two problems are given by 



X-dg,y{q')^K(G{0,q')-G{0) 



(12.2) 



X-d,,T{q')^>,TiO,q') 



(12.3) 



respectively. Solving (12.2) will follow as before but now we must make use of the diophantine condition on 
A G r^r rather than just on A. As before we set 



This leads to the conditions 



or equivalently 



C ■i + KEiQAX = KB-C -dq'X, A-C = kA 



C X^ 
X 



kEi(a 



kB-C ■ dq>X 



iA 



(12.4) 



Since we are interested in in solving the above equation for the unknowns f and Qai we would like the matrix 
on the left hand side to be nonsingular. Using the fact the matrix I is nonsingular together with a lemma 
comparable to le mma (5^) it follows easily the matrix on the left hand side of ( 12.4 ) is nonsingular. The 
solution of ([12. 3|) follows exactly as shown previously in the paper. Next we look at the solution of 



A • dq,Y = p, 



(12.5) 



where as before (3 is defined as 



{KB(q') - C{q') ■ e - C{q') ■ dq,X{q')) - {kB - C ■ ^ - C ■ dq,X) 



Since the definition of H^{p',q') has a quasi-periodic time dependence not present before, Bj{q'),C{q') and 
X{q') in terms of will change. Consequently the terms that make up /3 will change. For instance 



^{B,{q') 



B. 



E 



Furthermore, when calculating {C{q') ■ dq'X{q'))j three more terms must be added to this expression 



Previously these terms would dropout of the expression since dq' ^_^y{q') = 0. The solution of ( 12.5 ) splits 
into a quasi-periodic part and an exponential-with -respect-to-time part given by the following 



(12.6) 
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A • d,.T,{q') = Pf{q') = ^ PfAt)e''-'' G Ap^25,a-25, (12.7) 
where pf{q') and pf{q') are obtained from the expressions for KB{q') — kB, C{q') ■ ^ — C ■ ^, and C{q') ■ 



dqiX{q')) — C • dqiX. As done previously in the paper one can solve ( 12. € ) and ( 12.7 ) to obtain Yj{q') with 
all the necessary bounds. The rest of the bounds and the iterative process are obtained as before. The main 
result is described in the following theorem. 

Theorem 12.1 

Consider the Hamiltonian H{p',q') = H°{p) + KH^{p,q,t) + t = H°{p') + KH^{p',q') of 
{Ci,C2,Ci,C2,v.ji)p' ,q' -exponential form defined on B x T" where H'^{p') = H^{p) + r and n << 1 is a 
small perturbation parameter. Assume the perturbation has the following form 

H\p\q')^Q{p',q') + E{p\q'), 

Q{p',q')^ J2 sl{p'y''-^''''\ E{p',q')=J2hlip')el{t)e^^-'^. 

feeZ"+'" kei" 

Fix PqCzB and denote 

Assume there exists positive numbers F, 7, p, cr, d, k, all less than one and L, such that p > a and 

Aef^r, |A|<L, H°,H^ e Ap,a.,po, d\\i\\ < \\C*v\\ < d-^\\il Vi) e C", 

moreover, \\H\\p^a.po < 1> (the last condition will imply the isoenergetic condition on I). 

Then there exists positive numbers E,k' ,E' ,C,' , p' and a' , such that if \\H^\\p,a,pa < E one can construct 
a canonical analytic change of variables 

{V',Q') ^ <l>{V',Q') = {p',q'), 

with (f> £ Api ^lyi , which brings the Hamiltonian H into the form H' = H o (f> given by 

H'{V', Q') - {H o cp){V', Q') - a' + (1 + n'E'C)\ ■ V + 0{\\V'f), 

where a',(' G M. The change of variables is near identity in the sense that — identity\\pr^^ — > 
as ||-ff^||p,CT,po ^ 0- ^'^ '^^'^ coordinates Hamilton's Equations are given by 

Q' - ^(7", Q') = (1 + n'E'OX + OiWV'W), V' = Q') = 0{\\V'r). 

The solutions are given by 

Q'{t) = (1 + K'E'OXt + Q'o V'{t) ^V^^ {Vo, To) = 0. 

We can write out explicitly these solutions Q'{t) = (Qi(i), Qnit),T{t)) = (1 + n'E'OXt + Q'^ = {1 + 
K'i?'C')(A, l)t + Q'q which indicates in the new coordinate system the phase space T" x M x R" x R is foliated 
by invariant infinite cylinders each sustaining quasi-periodic motions identified by the frequency X{Vq) and 
evolving in the t direction. Note under the transformation, the time variable T{t) = (1 + n' E' C,')t has shifted 
by a constant proportional to the size of the perturbation. 
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13 Conclusions 



In this paper we looked at time aperiodic perturbations of nearly integrable Hamiltonian systems. To prove 

a KAM-type theorem we used the Lie scries formalism in an attempt to construct a scries of generating 
Hamiltonians from which we found a series of near identity canonical transformations. This series of canonical 
transformations is used to put the original Hamiltonian in a normal form which has a smaller perturbation 
at each step of the iterative process. Since we considered an aperiodic time perturbation, the Fourier series 
methods normally used to solve for the generating function no longer apply. We therefore made use of Fourier 
transform techniques which allowed us to understand how the analyticity domain of the Hamiltonian changes 
through the iterative process and how fast the perturbation shrinks. We considered the general form of the 
Hamiltonian H{p, r, q, t) = H°{p, r) + kH^{p, t, q, t) where r is conjugate to t. In the first case considered the 
perturbation decayed exponentially in time to a term depending only on angles. In the second case we added 
to this exponentially decaying term a time quasiperiodic term. In both cases we saw the time apcriodicity 
had to decay exponentially in order to have control over the analyticity of the Hamiltonian. Moreover, since 
the p', (/'-exponential form of the perturbation is essential in solving the partial differential equations that give 
the generating function, it was necessary to prove the p' , g'-exponential form of the perturbation is preserved 
after each step of the iterative process. 

A Properties of Analytic Functions on T"^ x C 

As mentioned before, techniques from multi-variable complex analysis are used in this paper. For exam- 
ple, a generating Hamiltonian with converging estimates can be obtained via Fourier methods assuming 
appropriate domains of analyticity for the perturbation of the Hamiltonian. Consequently wc derive some 
important properties of analytic functions on T" x C. We recall a common Fourier expansion for a complex 
analytic function depending periodically on the the variables q\,...,q2- Consider a multi-periodic function / 
analytic on some open strip of C" with periods Li, ...,Ln > 0, it is possible to represent / by the formula 
= X^fcez" c^e^'^'^i'^i'^"' where the coefficients Ck € C have the relation = c-k, and 

Cfc = lo' l^r-Zo"^" ^e-'(-^^^9^+-+-"'=''«")/(gi, ...,g„), where iv, = 27r/L„ i = 1, ...,n. 

A function also depending aperiodically on a variable t G C can be expressed similarly in a Fourier 
expansion. Let F{q,t) E Ap^r,,pa be real value, continuous, and 27r periodic in each qi. One writes F{q,t) = 
Ylik&^ fk{i)e"^'^, where fk{t) = j^Ty^ lo'' ' ' ' lo'' t)e~^'''^dqi ■ ■ ■ dqn- The following lemma gives a useful 
bound for the Fourier coefficients of a function F{q, t) depending 27r periodically on q and aperiodically on t. 

Lemma A.l 

For F e Ap^a and given F{q,t) = XlfeeZ" fk{t)e^''''^, then for every fc e Z" and for all t G C we have 
\fk{t)\ < ||F(t)||^e-l*=l^ where \k\ = l^il- 

Proof : From the previous argument we can write 

Mt) = I n<lu ...,qn,t)e-'^^ "^"^dq, ■ ■ ■ dqn. 

Since F S Ap^^^, we can shift integration paths as follows 

fk{t) = Fiq, - i^^p, g„ - ^^p, t){fl er''^^''-'^'''}dq,...dqn. 

Note in case ki = we don't shift the path in that direction. 
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Then the following estimate its clear 

m)\ < ||F(t)||,e-EI'=.l)P^^^e-'=-'d(?<||i^(i)||pe-l^-'''. 

(Note \\F{t)\\p is a function oft). □ 

The following lemma is useful to construct a function in the form of a Fourier scries with a specified 
strip of analyticity given the Fourier coefficients satisfy an exponentially decaying bound. It will become 
evident later how this will enable us to construct a generating function in the form of a Fourier series whose 
coefficients have a specified form depending directly on the perturbation of the Hamiltonian. 

Lemma A. 2 

Suppose for a positive real constant p <1 and every k G one has 

|/fc(t)| < C(^)e-l'=l^ C{t)>0 ,vteC, 
with fk{t) € Aa, supjj^{C{t)) = C, and consider the function 

F{q,t)= /feWe*'- 

feez" 

For any positive 6 with S < p one has F G Ap-s,a- and \\F\\p-s^a- < C (|)"- 
Proof : Let ||/m g|| < p — (5 then 



\F\ 



p—S,a 



sup 



feez" 



\k\S 



fcez" 



i! e z" \fe=0 / \ ^ / 



e z" 
fcl > o...fe„ > 



For any positive 5 < 1, we have < 2 which implies < |, from which follows that 



\\F\\p-5,a < C2" 



= C 



The bound holds for any finite sum of the terms fk{t)e^^''^ G Ap-s^a, o-nd any such finite sum is analytic 
in Dp-s,a- Therefore, F is the limit of a uniformly convergent sequence of analytic functions and thus also 
analytic. □ 

Lemma A. 3 

For any K,s,S > one has the inequality < (^)* e^''. 
Proof: For any a; e M we have x < e^~^ . Let x = then 



so that 



K < -^e^* or K 
oe 



■<-ii) 



66 



B Bounds on the Kolmogorov Normal Form in terms of the orig- 
inal Hamiltonian 



Proof of Lemma |5 

Clearly 



Assuming d\\v\\ < IIC*^!! < d \\v\\ where 



C* 



C* 




it follows 



|C*u|| = max 



Also, given v = {v, w) — {vi, w„, w) for v e C", w e C, and v e C"+-^ it follows 



and 



lC*u|| = max 



This implies ||C*u|| = \\C*i\\ and therefore ||C*u|| = \\C*i\\ < d-^\\d\\, Vv e C"+i and e C". 
Finally, since ||u|| = max(|ui|. |-D„|) < max(|ui|, |-D„|, |w|) = it follows IjC*!'!! < . □ 

We summarize some results and give some new definitions. We defined previously the (n x n) matrix 
C* = d^H"{0) and assumed d\\v\\ < ||C'*w|| < d-^\\v\\ i £ C". We defined the (n + 1) x (n + 1) matrix 
C* = dl,H°{Q) and obtained the upper bound ||C*u|| < d-^v\\ v e C"+i. We defined the (n + 1) x (n + 1) 
matrix C — dp,H°{0) + Kdp,H^{0, q'). We define the (n + 2 x ?i + 2) isoenergetic matrix 



1 = 

where A = (A, 1), A G C"+i and A G 
Proof of Lemma W- 



C* A^ 
A 



C* (A)^ 
1 



A 1 







1. max(P||p,.,||i3||p,.)<2f. 



We assume for the original Hamiltonian ||i?^||p,o. < E. Recall 



A{q') = (iJi(0,g') - B{q') = 

It follows ll^llp.cr < 2|l_ff-'^|lp^cr- Applying Cauchy's inequality and the fact p > a we obtain the bound 



67 



Therefore 



ll^llp,, < 2E<— (a<i), < :^ < ^ 



from which foUows 



2E 

max(P||p,^, < — =Ei. 

a 



2. \\Cv\\p^„ <iii-^\\v\\ 



By definition we have 



Therefore 



and it follows 



pp. TT p)2 ttO 



\\{C-C*)vh.„<{n + l)\\C-C*\\,Av\\ 



Applying Cauchy's inequality we obtain 



Recall from the statement of theorem 3.1 



kE : 



212A2 



One can verify (using d < 1 



from which follows 



±i2 f^V 

fT2 V32/ 



Now recaU ||C*i;|| < d-^\\v\\, Vi; G C"+i. Writing C = C* + (C- C*), we obtain 

IICHIp,. < ||c*HI + ll(c-c*>||p,.<(d-i + ^)||HI<2d-i||«||, 

taking m = |, the result follows. □ 

The next Lemma will be used in the proof of the Lemma on the Isoenergetic Nondegeneracy Condition. 



Lemma B.l Given two positive functions f{x) and g{x) then 

Tiia,yi{ f — g) > max(/) — max(g). 
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Proof: 

Let go be fixed. Then clearly 



Next set go = ma.x{g). It follows 
Consequently 
□ 

Proof of Lemma |5 



max(/ - go) = max(/) - go- 

X X 

go>9 and f ~ go < f - g. 
max(/ - g)> max(/ - 50) = niax(/) - 50 = max(/) - max(5). 



Let V = {v,vi,V2) where v € C^,vi,V2 E C. The following estimates follow 

= max (^\v2\, \X-v + vi\, | ^QjWj + V2Xi\j 



C* 
1 
A 1 



V 
V2 



- (1 H '^Ij^^ + W2 Ai 1) > max | (| ^ C-jVj \ - \v2X,\^ | > max (| ^ C-^v^ \j - max (\v2~Xi\j 
j 3 3 



> 



and we set I = ^\d — L\. □ 
Proof of Lemma |5. 
Recall 



Then clearly 



I = 



C A^ 
A 



C* A^ 
A 



> 



C{q') A^ 

A 



C* A^ 
A 



-Kdl,H^{Q,q') 




ndl,m{t)) \ 
/ 



> 



> 



A\v\\p,<y - 



id^p,H^O) 




P,<J 



> 



P,I7 







> 


kE 
9" 


) 


M\p,<y 




p,a 







where we use in the last line Cauchy's inequality and the fact ||77^(0, ^OIIp.ct < -B. □ 

C Lemmas for Canonical Transformations 

In what follows we present several estimates relevant to canonical changes of variables defined by the Lie 
series method. For x{p',q') S ■^p,cr define 



Ap,( 





dx 




dx 




( 


dp' 


p,cr 


dq' 
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Lemma C.l 

For any x, / € Ap^^ o,nd positive 6 < a < p we have 



1- \\{X,f}\\p-S,a-5<2{n+l)(^^^ 

2- \\{x,{x,f}}\\p-s,a-6<8{n+l) 



2 / ^p,o- 



Proof : 



||{x,/}||,-5,.-.<2(n+l)(4^) 



By definition 
It follows 

\{xj}\ < E 

i 

which implies 



{x,/} = E 



9x 9/ _ dx df 
dp', dd dd dp'. 



dx 



df 



dq'i 



+ 



dx 



df 



dp'i 



< 



dl 



||{x,/}||p,.<2(n + l)(^)l|/|U. 



df 



<x;/-^\\f\u. 



\\{x,{x,f}}\\p-s,a-5<8{n + ir(^^y 



,\P,<7 



We have by definition 



{x,{x,/}} = E 



dp'.dq'^^'^^ dq'^dp'^^'-' \ 



= E 



dx d 
dp'i dql Y 



dx_dl_ _ dx df 
dp', dq'j dq'j dp'j 



dx d 
dq'i dp'i ^ 



d]^df_ _ dx df 
dp', dq'. dq'. dp'j 



E 



dx 
dp' 



d\ df ^ dx d^f d\ df dx d^f 
dq'idp'. dq'j dp'j dq'idq'j dq'idq'j dp'j dq'^ dq'^dp'- 



dx 
'dql 



d\ df_ dx d^f d\ df dx d^f 



dp'idpf^ dq'j dp'j dp'^dq'^ dp'jdq'j dp'j dql dp'^dp'j 



There are eight expressions in the sum, each can be estimated using Cauchy's inequality. The first two 
expressions are estimated below and the rest can be estimated the same way. 



dx d'^x df 



dp'j dq'jdp'j dq'j 



< 



p—S,a—5 



x: 



p,a 
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dp'i dp'j dqldq'^ 



< 



Xn,(T 



p—S,(T—5 

Each of the eight expressions in the summation represents {n + 1)^ terms. Therefore the estimate gives 

,'y* ^ ^ 

\2 / ^p,"' 



\\{X,{X,fmp-5,a-5<^{n+lY 



□ 



We now review the Lie series method in which the desired canonical transformation (p is given by the time 
one flow of a Hamiltonian x called the generating function. If one evaluates a function / on the solutions 
of Hamilton's equations generated by x o'^^ can write ^ = {X;/}- That same functions / defined on 
the phase space is transformed f ^ Uf = f o (p. Therefore this transformation of / can be written as 
Uf = f + {x, /} + ^{Xj {X)/}} + ••• which is exactly the Taylor expansion of /. The transformation of a 
function can thus be expressed entirely in terms of the generating Hamiltonian x without knowing the flow. 
The following lemma gives the conditions under which such a canonical transformation exists and several 
useful estimates. These estimates are key in the iterative lemma where the size of the perturbation after 
applying the canonical transformation must be known. 

Lemma C.2 

Suppose x(p')?') € Ap^a- o-nd its derivatives also belong to Ap^a for given positive numbers p,a. Consider the 
corresponding Hamiltonian system 



^ dp' 



Assume that for some positive 5 < a, 5 < p 



Xp^a 



max 



dx 




dx 


dp' 


p,(T 


dq' 



s 

<2- 



Then, for all initial conditions {P',Q') € Dp-^^^-s, the image of a point {P',Q') € Dp-s^^-^ is contained in 
Dp^a , thus defining a canonical analytic transformation 



The operator 



■ Dp-S,a-S Dp^„, 4> G Ap-S,a-S- 



f^Uf = focP, 



is then well defined and one has the estimates 



1. identity \\p-s,a-5 < Xp,^' 

2. \\Uf\\p_s,a-5 < \\f\\p,a, 

3. \\Uf - f\\p-s,.-5 < 4(n + 1)4^11/11,,,, 

I \\Uf - / - {x,f}\\p-8,o-6 < 32(n+ 1)2 (4^)' 11/11,,,. 
Proof : 
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One first shows the existence of solutions through an arbitrary initial condition (P',Q') G Dp^s.a-S- 
This follows if the vector field is Lipschitz throughout Dp^s,i7-5- Using Cauchy's inequality in the domain 
Dp_s/2,a-s/27 the Lipschitz constant for the vector field 

dq' ' dp' ' 

can he estimated as 



dq'dp' 



p-S/2,(j-S/2 



< 4(71+1)^ 




Next we argue that this local solution exists for t ^ \, and does not leave Dp ,j in that time. By assumption, 
the maximum velocity of the vector field is hounded hy 

max(||p'||,||q'|l) < x*p,a < 2- 
Hence the maximum distance that a point in Dp^s^„-s can move in time t — 1 is ^. Therefore 



6* : Dp^s,cr-5 Dp_s/2,cr-5/2 C Dp^s, 



Next we prove the estimates. 



1. \\(f)- identity \\p-s,a-S < Xp^ 



This follows from the mean value theorem 



For 0<t<l, 110— identity \\p^s,a-s = \\ 



,t=l it=0 



||p-5,.-^<||<^*||(l-0)<x; 



2. WUfWp-S^^-S < \\f\\p..a. 



This is trivial, \\Uf\\p-s,a-s = \\f ° 4>\\p-s,a-s but (j) : Dp_s,a-5 — > Dp,s which implies \\f o(l)\\p_s^a-s < ||/|Ip,<t- 



3. \\Uf-f\\p-s^,-5<^{n+l)^\ 



By definition we have Uf — f o cj)* ^ . Taylor expanding in t gives 



from which follows 



As argued ahove. 



/o0* = /o0O + -(/o0*)|,=,,, 0<t'<l, 



Uf^f + f^U^t' or Uf-f = ^^\,=t,={x,f}\t=t' 



<^(F', Q') e Dp_s/2.a-S/2, 0<d<l, (P', Q') e Dp^s.aS 



Hence, using the estimate in Lemma C.l, with 5 in the Lemma replaced by ^, we obtain 



\\Uf^f\\p_S,a-5<^n+l) 
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4. \\Uf-.f-{xJ}\\p-s..-s<32{n+l)^^^) \\f\\p 



Taylor expanding in t 



and applying the mean value theorem gives 

Uf - f - '^\t=t" = 2 dZ^'*"*" ^ 2^^^' /}}!*=*"■ 

As above, 

<P{P', Q') e Dp_s/2,a-s/2, 0<d<l, (P', Q') e i?p_,-, 

Using this and Lemma C.l, with 5 in the Lemma replaced by |, gives 



\\Uf-f- {X, mp-5,.-8 < 32(n + 1)2 ( ^ V ll/llp, 



□ 



D P.D.E on a Torus 

Proof of Lemma |87 

We write the coefficients of the Fourier series as fohows 

g,^jl^Jl---iyig)e---^d,,...dg,.. 
Since G G Ap, we can shift the path of integration as follows 

(If fcj + we don't shift in that direction) and \gk\ < e^l'''l''||G'(q)||p. □ 



Proof of Lemma |8. 

We begin by expanding F and G in Fourier Series 



F{q) = /fe^'" '' ^(9) - E 9ke- 



ik-q 

k£Z" fceZ"\o 

substituting in ( ^.ll ) gives 



i fceZ" fceZ"\o 
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or 



which impUes 



fk - t|^, k e Z"\0. 
i\ ■ k 



X! 9ke 

fc6Z"\0 



ik-q 



From this expression we see why we must require G = in order to construct a solution. Recall that A G fir 
so that A • fc 7^ 0. We next prove convergence of 



9k ik; 



iX ■ k 



From lemma B.l since G G it follows \gk\ < ||G||pe~l'^l'', and the following estimates, using ||fc|| < |fc|, hold 

|A-fc| >r||fc|r", |A-fc| >r|/cr". 



Using lemma A. 3 gives 
Hence, we use these estimates to estimate the Fourier coefficients as follows 



\k\5 



\fk\ = 



9k 



X-k 



< 



\G\\ 



.|^.|ne-|fc|p < Mk ( M%-|fc|(p-5) ^ c-e-l 
r Veo/ 



k\{p~5) 



where 



C 



F \e6 



Using lemma A. 2 replacing p with p — 6 we obtain 

Next we get the first estimate on the norm of F in terms of the norm of G. The previous estimate can be 
rewritten as 



\F\\p-25 < F- 



1 f 16«\" ii^ii _ /167i\" \\G 



We have 



smce 



w ^2 



An+l 



n + 1 



|G||p = 



F(2<5) 



2n ■ 



n+1 



> 



1671 

e 



n+1 



ri+l 



Using this inequality, the estimate can be written as 



> 



Letting S ^ ^ gives the result. 
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iX ■ k ' 



Finally, we prove the last estimate in the lemma. Recall that 

fcGZ"\0 

then, formally, we have 

■' fceZ"\o fceZ"\o 

We denote the Fourier coefficients of ^ by hk^ = ^(fcj)i j — l,---,'n. Therefore l/i^^. | < 
r-i||G||p|fc|"+ie-l'=l''. From lemma ^ we have 



n + l 
eS 



Form lemma L4.2 we have 



where 



so 



dF 



dq 



< C 



P-2S 



n+l 



1 /n+l 



n+l 



IIGIlp, 



dF 



dq 



< 



4" /n+l 



p-25 



n+l 



|G|| 



Let (5 ^ |, then the estimate takes the form 



dF 



< 



p-5 



r(5"+i 



IIGIlp. □ 



E Fourier Transform Lemmas 

We present some important results concerning a function and its Fourier transform defined on complex 
domains. 

Lemma E.l 



Suppose f{z), z = X + iy, is analytic for —a < y < b where a > 0, 6 > 0. 
In any strip in the interior of —a < y < b, let 

O {e-^^^'^y-') {x ~* oo) 

O (e(^-^)^) {x -oo) 
for every positive e, where v > 0, fi > 0. Then there exists a 5 > e such that F(u)), oj — u + iv, defined by 



F(c.) = -= / /(C)e'^-^dC, 



is analytic in the strip — /i + 5<v<i'~5, satisfies 

0(e-(^-^'") (u->oo) 



F(u;) 



O (e^-^-^)") {u -oo) 
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and 

fiz) = -^J°^^F{w)e"'^du;, (E.l) 

for every z in the strip —a < y <b. 
Proof: We have 

to find the domain of analyticity of F{uj) we consider the sequence of functions {-F„(a;)}^Jj" where 

1 f" 

Fn{co) = -^j^J{Oe-'<"dC 



and show that converges to -F(w) uniformly. Consequently, since each of Fn{u)) is analytic in 

the strip —ji < v < v, uniform convergence will imply F{uj) is analytic in some strip to he determined. Let 
n = inax(ci, |c2|) where ci and C2 are the constants in the definition of a function of exponential order in the 
real limit. Then, taking integrals along the real line. 



\F4co) - F{co)\ 



< 



i/(c)iie-<-idc+ / i/(c)iie-*^n'^c 

J — oo 

1/(01 |e-^<"||e^''|dC+ r i/(C)l|e-^^"iie^''|rfC 

J —OO 

POO pn 

I |/(C)|e«''dC+ / |/(C)|e«''rfC 

n J — oo 



'2tt 
1 

1 



< 



2tt 
1 



27r 
1 



poo pn 
J n J —oo 



-{v — e — v) 



+ 



^ (^,-e+v)C 



jJi — £ + V 



The maximum finite hound for the expression above can he obtained by picking for the zeta intervals [n, oo) 
and (—00, — n] , v = v — 5 and v = —ji + 5 respectively where 5 > e. Note that, if v is picked anywhere between 
{iy' — 5, v\ or [— /i, — /x + (5), the expression above becomes unbounded. On the other hand ifv is chosen anywhere 
on—^, + 5<v<y — 5 the expression above is bounded. Therefore 



\F^{u:)-F{u:)\< 



2C3 e-f''-^)'^ 



< e 



(^-£) 

for —fi + 6<v<i' — 6 where C3 = niax(Ci, C2). We can solve for n 



-{S — e)n < In 



2C3 



{5-e)e 



so 



n > In 



27r 
2C'z 



i6-e)e 
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We may take 



N{e) = In 



2C. 



{S - e)e 



Therefore, for any e > we have defined a iV(e); independent of ixi, such that \F{lo) — F„(w)| < e, for all 
n > N. Consequently the series {Fn(u>)} converges uniformly to F{uj) for —fj, + 6<v<i' — S. Hence F{u)) 
is analytic in this strip. Next we prove the order results for F{u)). By applying Cauchy's theorem we may 
take the integral along any line of the strip parallel to the real axis. Thus 



\n^)\ 



< 



-1= / /(e-«??)e-*(«-^'')(''+*'')(i^ 

/2'K J-oo 

-7]U poo 

/ \f{^-iv)\e^''d^ 

J — oo 



'^"||ef''||e-''"||e-*''''|d^ 



2lT 



-r}u 



2-K 





|/(C-z7?)|e«X 



\f{i-i^)\e^''di 



+ / l/(^-ir?)|e«''dC+ r |/(^-ir?)|efdC 

C2 J —OO 

POO PC2 

(v) +Ci e-^-'-^^^e^^rfC + Cc2 (v) + C2 / e^^'-^^^e^^'d^ 

t/ ci —00 



277 



-rju 



-{p — e — v) 



00 (726^''"'^+"^^ 



ci iJ. — e + V 



where 



Ccriv) = r |/(^ - ir?)|e«''dC, C,,{v) = H - iv)\e^''d^. 
Jo Jo 



For the interval [ci, 00) pick v = v — 6 with e < 5 such that 



-{v-e-v)(. 



-{v — e — v) 



< 



-(5-e)ci 

5-e ' 



Similarly for the interval (—00, C2] pick v = —fj, + 6 and with e < 5 we have 



oC/n-e+w)^ 



fj, — e + v 



< 



,{5-e)c2 

5-e ' 



We bound Cci{v) and Cc2{v) by choosing v = u — 5 and v = — /x + 6 respectively so that 



Jo 

CM < I \f{^-iv)\e^-^+'^^d^ = C,,. 

J C2 

Consequently, with cs = niin(ci, \c2\), C3 = max(Ci, C2) and Cc^ = max{Cci , Ccj) 



\F{u;)\ < 



2n 



(<5-e) 
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and by taking rj arbitrarily near to —a and b the order results for F{lo) follow. (E.l) can be proved directly 
by the theorem of residues. Let —a <a<y<f3<b. Then 



— r 

\/2^ Jo 



Similarly, 



F{Lo)e''''^dLo 



1 

1 

27ri 



1 

2^ 



--ia+oo 



/(C)dC / e-^^^^-^'^'^du 

Jo 



/(C) 



1 1 f(c) 



lil3-oc ^ ~ C 

Now we consider the following counter clockwise contour integral along a rectangle 



/(C) 

c C - ^ 



/(C) 

i/3+p C ~ ^ 



dC= / ^dc- 



/(C) 

c-^ 



/(C) 

c-^ 



ia.—p C 



/(C) 
L2 C- 



dC+ / :f^dC 



/(C) 

c-^ 



a Li 



We show the contribution from each of L2, L3, L^, Lq is zero. For L2 we use the following parameteriza- 
tion C(y) — p + iy; —a < y <0 and since f{C,) must be bounded for all ^ by some real number A 



/(C) 



so 



lim 

p — >QO 



/(C) 



< 



A 



\p + iy- z\ 



dc 



AoL 



< lim 

P^oo \p + iy — z\ 



A similar procedure with Q <y < (3 in the parameterization gives 



lim 



/(C) 

c-^ 



dC 



< lim 



Aa 



p^oo \p + iy — z\ 



= 0. 
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For L5 we use the following parameterization C,(y) = —p + iy, < y < /3- 



/(C) 



so 



lim 



no 



< 



A 



\-p + iy- z\ 



Aa 



< lim T 

P^oo \- p + iy- z\ 



= 0. 



A similar procedure with —a<y<0 in the parameterization gives 



Therefore 



and 



□ 



lim 



/(C) 



dC 



Aa 

< lim : - 

p^oo I — p + ly — z\ 



= 0. 



JC Q - Z Jij3+p Q - Z J-ia-p Q - Z 



f{z) 



F{uj)e-"'^dijj = 

27r./-^ 27^^ 



i/3+oo 



i(3 — 00 



/(C) 
z-C 



dC 



/(C) 

C-z 



dC 



1 

2^ 



/(C) 

c-^ 



— [27ri/(^)] = /(^). 



Lemma E.2 

Suppose F{oj), oj = u + iv, is analytic, regular for —u < v < p, where 
v > 0, p > 0. In any strip in the interior to —u < v < p, let 



F{u;) = 



O (e^"-"^)") {u -00) 
for every positive e, where a > 0,b > 0. Then there exists a 6 > such that f{z), z = x + iy, defined by 



1 f°° 



is analytic for —b + S<y<a — 5, satisfies 



O (e-(''-^)^) {x 00) 
O (e^''-'^)^) {x -00) 

F{uj) = J^J{z)e-'^^dz, 
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for every oj in the strip ~v < v < jj,. 



Proof: We have 



1 r°° 

f{z) = -j= / F{Oe^^<dC, 



and define the sequence {fn{z)} where 

fn{z) 



F{Oe^<dC. 



Similar to Theorem E.l we can show that this sequence converges uniformly to f{z) for —h + S<y<a — S. 
And since each of fniz) is analytic for —b + 6<y<a — 6, it follows that f{z) is analytic in this strip. Next 
we prove the order results on f{z). Applying Cauchy's theorem we may take the integral along any line of 
the strip parallel to the real axis. Thus 



l/WI 



< 



< 



1 



— OO 

oo 



V2tt 

,-xri 



1 f°° 

|F(C + i77)||e"«||e-"''||e-^«||e-'^''|de 
\F{i + iTi)\e-y^di 

POO 

\F{^ + tij)\e-y^d^+ / \F{^ + iTj)\e-y^d^ 



/27r 

— XT] poo 



+ / \F{^ + iTi)\e-y^d^+ / \F{^ + iT^)\e~y^di 



-xrj 



< 



/2Tr 

— XT] 



e-(ft-e)«g-y? + Cc, (y) + C2 / e("-^)«e-^«dC 



-(b-e + y) ci (a-e-y) 
To bound Cci (y) and (y) we chose y = —b + 5 and y = a ~ 5 respectively so that 

Jo 



For the ^ intervals [ci, 00) and {—00, C2] we pick y = —6 + S and y ^ a~ S respectively with S > e and obtain 
with C3 — min(ci, \c2\), C3 = max(Ci,C2), and Cc^ = max(Cci , Cca) 



iS-e) 



and by choosing rj arbitrarily close to —v and /i the order results follow. The rest of the theorem can be 
proven with the theory of residues exactly like the previous theorem except that the contour integral will be 
taken clockwise. □ 
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Lemma E.3 

Given G{q, t) G Ap with 



then 



Proof: 



9k{t) 



9k{t) = 



G{q,t)= ^ 9k{t)e'''' 



e)(e(«-<^)t«) (tR -oo) 



0(e-(f'-e)tH)0(e-|fe|p) {tn ^ oo) 

e)(e(«-^)*H)0(e-l*=lP) (f^ ^ -oo) 



(27r)^ 



TT •/ — TT 



G{q,t)e-''''idq 



where the integral is along real axis. Lift the integral as follows 



9k{t) 



1 



TT /•TT 



„-|fc|p /■TT /■TT jr. jl, 

5fe(i) = J ■ j G{q, - ij^^p, g„ - i^p) • e-'^ Uq 



\hri 



and \gk{t)\ < e-\''\P\\G(q,t)\\p. We know 



\\G{q,t)\\p= sup \G{q,t)\. 

Let q* be arbitrary in the interior of the strip \qi\ < p. Then by the above and the assumption on the order 
ofgk{t) 



\G{q\t)\ 



E 9k{t)e^'''' 



feez" 



\k\ip-q'j) 



feez" 



feez" 



feez" 



fcl , ...kn > 

= C'e-^*«2" I^E'^^^'^'j =Ce-^*«2" ( i_"^g-e ) ' 

where A is b — e or —a + £. □ 

The following lemma shows that under certain conditions we are assured the existence of the complex 
Fourier transform. 
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Lemma E.4 

Suppose f{t), t = tii + iti, is analytic for —a<ti<b where a > 0, 6 > 0. In any strip in the interior of 
—a <ti<b, let 



f{t) = 



O (e-(''-^)*«) {tR oo) 



O (e(''-")*«) {tR -oo) 
for every positive e, where v > 0, jj, > 0. Assume the Fourier transform exists 

1 r°° 

Then the complex Fourier transform defined as 



< oo. 



F%u;) = -= / f{C)e-'^^dC, 



with —a < <b exists and F'^{u) = F{uj). 

Proof: We consider the counter clockwise contour integral on the rectangle R show below 







tl 








L2 












-p 




P 



tR 



Since f{t) is analytic for —a < tj < b we know the contour integral is zero. Consequently 

-I -hhl ^ 

We analyze the integral along L3 using the parameterization t = p + iti such that 

[ fk{t)e-'^*dt = f fk{p + iti)e-'^^''+'"hdti= ['^ fk{p + iti)e-'^'e^Pe^*'e'^*'idti 
J L3 Jo Jo 



and 



Note although 



[ fk(t)e-'^'dt 




J L3 





\fk{t)\ < Cie-(^-^)*« ci<tR< 00, 
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we can choose a constant C3 large enough such that 
and 

< r C3e-'^''-'^Pe''Pe''''dti. 



[ fk{t)e-''^'dt 


<-L 


JL3 





By Theorem E.l, F{uj) is analytic inside the strip —ji + 5<v<v — 5 where e < 6 so that 

fk{t)e-"^*dt < C3e-'^''-''>Pe^''-^'>Pe^"^f^f3. 



L3 



Clearly 



lim 

p — >oo 



fkit)e~''^'dt 



L3 



= 0. 



Similarly for L4 using the parameterization t = —p + iti 

fk{t)e-'^'dt = f f,{-p + itj)e-'^^-P+'''hdtj^ r fk{-p + iti)e'"'e-^Pe^''e'^''idti 
and 



1 fk{t)e~'^'dt 


< j 


J Li 





With a large enough constant C4 



and clearly 



f h{t)e-'^'dt 


< j 


J Li 





lim 

p—*oo 



fkit)e-''^'dt 



Li 



= 0. 



It follows 



lim — / = lim 

P^°° Jl2 P^^JLi 



1 roo 1 pco+ip 

F{u) ^ -= / /(C)e--«dC - -= / fiCy-^dC = F^(c.). □ 

F Properties of p' ^ g'— Exponential Form functions 

In this section we study some characteristics of functions of the form 

F{p',q')^f{p') + r{p',q')+g{p',q'), 



ik-q 



g{p\q')^ J2 hk{p')ekit)e 



ik-q 
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where F{p' , q') G Ap,a and ek{t) is of exponential order with respeet to time. We will refer to these functions 
as functions of j?', g'— exponential form. The Hamiltonian under consideration will be of p', g'— exponential 
form. There are several things we want to know about these functions to carry out the iterative process 
of the KAM proof. First, we want to be able to extract from F(p',q') the quasiperiodic and exponential- 
order-with-respect-to-time parts as functions of F{p',q') itself. The importance of this will become evident 
in the next section. Also, we want to show the derivative with respect to p'j or q'j, j = 1, ...,n + 1, of a 
p' , (jt'— exponential- form function is another p' , g'— exponential-form function. First we prove an important 
characteristic of g{p',q'), the function which is of exponential order with respect to time. 

Lemma F.l 

Given a function g{p', q') G Ap^a of the form 

9{p',q')= Yl hk{p')eu{t)e'^-\ 

where ek{t) is of exponential order with respect to time, g{0) = 0. 
Proof: We simply show the time average ofg{p',q') is zero 

lim i r g{p',q')dt = V /ifcb') lim 1 C eu{t)dte''-^ = 0. □ 

Next we show how to express g{Q, q') in terms of F{p' , q'). 



Lemma F.2 



Given a function F{p' , q') e Ap^a of p' , q' -exponential form as defined previously, T{Q,q') = F{0,q') — F{0) 
represents the q- quasiperiodic and exponential-order-with-respect-to-time parts of F{0,q'). Moreover T{0,q') 
has zero average. 

Proof: First 

F(0) = 7(0) + f (0) + ^(0) = /(O) + so(0). 

Therefore 

F{0,q')-f{0)=f{0)+r{0,q')+g{0,q')-fi0)-so{0)= ^ Sfe(0)e^'=« + ^(0, g')- □ 

feezn\o 

Next we will show the derivative with respect to p'j or q'j, j = 1, ...,n + 1, of a p' — 5'— exponential- form 
function is another p' — g'— exponential- form function. It is easy to see this is true for the derivative with 
respect to p'j j = 1, n+1. 



Lemma F.3 

Let F{p',q') be defined as before. Then §^{p',q') j = l,...,n+ 1 are of p' ,q'- exponential- form and 
§^{p',q')&Ap,,. 

Proof: First we write out the derivative 

dF{p',q') ^ dm , dr{p',q') dg{p',q') ^ dfi^ . V ^"^(P') ^^-i . V ^^^^^'^ (fU^-'i 
dp', dp', dp', ^ dp'j dp'j dp'j +^4; dp'^ ■ 

Furthermore we know , as consequence of Cauchy's Integral formula, that if a function is analytic in some 
domain then all its derivatives are analytic in the same domain. □ 
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The next lemma show the derivative of a p' — g'-exponential-form function with respect to q'p j = 1, n, is 
another p' — g'-exponential-form function. 

Lemma F.4 

Let F{p',q') be defined as before. Then §^{p',q'), j = l,...,n, are of p' , q' — exponential-form and 
§lr{p',q')GA,,.. 

Proof: First we write out the derivative for j = 1, ...,n 

Furthermore we know . as consequence of Cauchy 's Integral formula, that if a function is analytic in some 
domain then all its derivatives are analytic in the same domain. □ 

Recall a function ek{t) is said to be of exponential order with respect to time if the following holds 

0(g-(i.-<r)tR) . < Ci < tij < oo 

C>(e(''-^)*«) ; -oo < tij < C2 < 0. 
The following lemma examines the derivative of a p' — g'— exponential- form function with respect to q'n+i — t. 
Lemma F.5 



e/c(i) 



Let F{p',q') be defined as before. Then for some positive S > a, ^{p',q') € Ap^a-S is of 
p' , q' — exponential-form. Moreover for the nth derivative we have the following estimates 



for all t e Vp „_s. 



Proof: We will prove this using Cauchy 's Integral formula. Since ek{t) is exponentially small with re- 
spect to time we have 

\ek{t)\ < Cie-^"-")*", < ci < ifl < oo, 
\ek{t)\ < C2e^^-'^'^, -oo<tR<C2< 0. 

Since ek{t) is analytic in the interior of the strip —a < tjm < a- we have 

where C is the circle of radius S centered at t. We used the parameterization w = t — 5e^^ < 6 < 2w and 
obtain 

et\t) = ^[ r^^^d^= r '''('^-^''f + '^^;-^''^^^h se^'d0 
^' 2-Ki jc - io ((5e»^)"+i 
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and for < ci < tR < oo 



Therefore 



/or a/i t e Vp^cr-s- Similarly for — oo < fij < C2 < with the parameterization u = t + 5e*^ we obtain 
□ 

We want to show the product of two quasiperiodic functions is quasiperiodic. 
Lemma F.6 

Let F{q), G{q) £ Ap be periodic with respect to q GR. It follows F{q)G{q) = H{q) is periodic with respect 
to q and H{q) can be expressed as 



Jkq 

feez 



where Ck is defined in terms of the series expressions for F and G. 

Proof: First we prove that the series expression for F and G converge absolutely to the functions for 
q GVp. Given 

Fiq) ~ ""ke'"' e ~ E ^'^e'"' ^ 

feez feez 



we have 



\ak\ < ||F|Le-l'=l^ 



Furthermore, given \qi\ < p — 5 for some positive 5 
Therefore 

|afce*'=«| < ||i^||pe-l^l''el'=l(''-*) < ||i^||pe-l^l^. 

A similar estimate holds for \bke^'^'^\. Consequently, '^^ez '^kS^'"^ and X^^gz bke^'^'^ converge absolutely to F{q) 
and G{q) respectively for q GVp. Next we look at the product H{q) given by 

/ \ / \ / k=oo \ / k=oo \ n=oo / k=oo \ 



Kk=—oo / \k= — oo / n=—oo \k= — oo 



Vfeez / Vfeez / \fe=-oo / \ 

= E ( E ) e"' = E C"^"" ^^P- ° 

nez Vfcez / nez 
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A similar result holds for quasiperiodic functions and can be proven similarly as above. Given 

F{q) = J2 ^ke'''-^ G{q) = ^ &fee^'=-«. 



It follows 



H{q) = F{q)G{q) = ^ ( E «"6,e» « ) e^'=" = ^ Cke'"-'. 



Lemma F.7 

Let F{p', q') e Ap^a be of the form 



F{p',q')=f{p')+r{p',q')+g{p\q'), 

r{p',q')= Mp'V'"'' 
9{p',q)= E hk{p')ek{t)e^'', 



where ek{t) is of exponential order with respect to time. Let xip'j l') G ^p',<t'; p' < P O'^d a' < a, be of the 
form 



x{jp',q')=X{q')+^-q' + Y{q')-p', 



where 



X{q, t) = y{q) + T{q'), y{q) = ^ y,e^'=-«, T(g') = ^ Xk{t)e"' '' 



fcez" 



feez" 



w/iere ^ is a constant vector. It follows, given some positive d < a', {X) € A.p'-s,a'-5 is of 
p', q' -exponential form. 

Proof: 

We write out the expression 



n+l 



{X,F} = J2 



9x dx dF 

dp'j dq'j dq'j dp'j 



n+l 

E 



dr 

M 



dq', 



dq: 



l7(^')+E|^(^')e--+E|7(^^')e.We-- 



n+l 

E 
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n+l 

E 



dq' 



(^(.').E||.Oe-.E|j.O.(.e-); 



n+l 

E 



5,•(g)^(p^gO +.^.(9')#7(y,9') 



dq' 



'dq', 



( dy ( .df , dy ^ dsk . 



dq[ 



3 fcez' 



dq'. dp'. dq'. ^^dv', dq', dp', 



VKe now mwst wnie out each of the nineteen expressions above and rewrite the entire expression in 
p' — q'— exponential form. The first term 



n+l „ n+l 

i:^M§^iP',<i') = E 



Vfcez" J 



E 



.fcez" 



feez" 



= E 



.feez" \mez" 



E I E E^'5-.^^^-^'^(^^')e^'"'p^''' 



E E E^'^'"'j'^j'*"'-™(^') I 

u;eZ" \m6Z"i=l 



iw-q 



The second term 



n + l n+l 

Y.sM)^iP',<i') = E 



9g 



.feez" 



J Vfeez" 



= E 



E Skje'^-'i E ikjhk{p')ekit)e"'-'' 



feez" feez" 
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E 



5^ iSm,ikjhk{p')ekit)e'^-'' e^'=" 



dt 



feez" VmeZ" 



ujGZ" \meZ''j=l 



T/ie third term 



E 

3 = 1 L 



E 



E E E^^'».iW^^-^'^(?'V'"M«*'= E E E*^-.^w("^-"^)i«— (^^ole^"-'. 

feez" \ mez" j=i / toez" \mez"j=i 



T/ie fourth term 



3=1 



1 \feez" 



= E 

n 

= E 



^kAty^'" E ikjhk{p')ekit)e"'-' 



.kei 



E E i^mjmjhkip'jekitje'^-" e^'^" 
.feez" VmeZ" / 



feeZ" VmeZ" 



E I E Y'^rn,3mhk{p'y{t)e'"''' I e^'=-«+ X ( X J^m,n+i{t)hk{p )^{t)e"^ A e*" 



feeZ" \ meZ" j=i 



feeZ" \mez" 



r/te ^/i/t term 

n+l 



E 



n+l 



d 



E 



.feeZ" 



T/ie sixth term 

n+l 



E 

i=i 



^Pi 



n+l 

E 

i=i 



#7 ( E y^^'^A E l^b')e^'=-' 



89 



= E 



E 



E ( E -.^"^fe(^^')e"" ' ) ' 
feezi Vmezi J 



T/ie seventh term 

n+l 



dp', 



E 



ri+l 

E 

i=i L"'^J 



1^ E.-*" ,i:SWe.We*. 



feez" feez" 



E 



.fcez" WeZ" 



fcGZ" \ mGZ" j=l 

The eight term 

n+l 

E 



= E I E E^-iy-|r(p')e.(i)e^'-« I e^'-'= E | E E *-i2/™^(p')e.-™(i) | e^""' 



n+l 

E 



1^ E ^f^') 



n+l 

E 



= E 



E ^kixu{t)e^'^■'^^{p') 



n+l 

E 



^7 E-^w^*ME|7(^')e-« 



= E 



E ^w^"-^ E £^^p') 



n+l 



= E 



feeZ" \mez" 



E I E E^-.-"^wS(^^')e^'"-|e--+E ( E ^^'^^yy""']''^ 



fcez" Wez" 



= E I E E^-.--W%^(P') E ( E %'(*)^(^'') 1 e^'"-"- 



t«eZ" \m6Z"i=l 

n+l 



E 



n+l 

= E 



5,^ 
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The eleventh term 



The twelve-th term 



= E 



J2 ik,x,{t)e^^-<^ J2 »^(p'>>^ity'' 



E 



.feez" WeZ" 



feez" \mez"j=i 



+ E f E ^w^(^>^w^^'"") ^"-^ 

= E I E E ^^j^rnit) ^^Q~r {P')ew-m{t) | 

+ E f E 



n+l 

E 



n+l 

E 



E |Ee.||(p')|e-- 



r/ie thirteenth term 



n+l 

E 

r/ie fourteenth term 



' n+l 



E E0 3-fb')e.W|e^'=-. 



n+l 

E 



n+l 

E 



E 



n+l 



E£. E^Me-0^;|^(.') 



n+l 



n n+l 

= E IEE%^MP;^(P') 

fceZ" \ j=i ;=i 



dp', 
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The fifthteenth term 



n+l 

E 



.1=1 ^0 fcez" 



n+l 

E 



E|.(E'5Me--)^';E|^(^')e-^ 



= E 

n 

= E 
= E 

feez" 

= E 



'n+l 



.;=i I 

n+l 



E ( E ( E -,^..pj||(p')e^-^ 1 



=1 VfceZ^ VmGZ' 
n n+l / 



EE( E «"^i'5'".'Pj|^(p')e' 



im-q 



j=i 1=1 Wez^ 

n n+l / 



EE( E i"^j^rn,ip f^zn p') 

j=l 1=1 VmeZ" 



The sixteenth term 



n+l 

E 



n+l 



E|7(^ME|^(^')e.We-^ 
;=i feez'i 



n+l 

E 



n+l 



E|r E'5^.'^"M^;El^(^')e^w^ 

;=i Vfeez" 



= E 



n+l 



.1=1 feez" 



feez 



= E 
= E 



n+l 



E ( E ( E ^rn,Srr.M^{p')ek{t)e-'' ] e 



ikq 



j=i li=i Vfeez^ \meT' 

n n+l 



EE( E i^i'Sm,ip'i^ip')ek{t)e' 



im-q 



= E 



j=i i=i 

n n+l 



EEI E imjS^^ip[—^^{p')e^-m{t) 



j=l 1=1 VrneZ" 



r/ie seventeenth term 



n+l 

E 



n+l 

E 



n+l 



n+l 



;=i 



E 

n+l ^ TT -P 
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= E 



fcez" 



n n+1 Qj 

j=i 1=1 ^-i 



Jk-q 



kei." \i=i "^^"+1 / 



The eighteenth term 



n+l 

E 



n+l 



.1=1 "3 



The nineteenth term 



n+l 

E 



El7fE-^MWe--).;E|7(^')e"- 



= E 



n+l 



n+l 



E E ^we*«.; E 



i=i feez" 

n fn+l 



feeZ" 



E 



n+l 



E E ( E ^m,^„,(t)pj||(p')e^'-^ 1 e^'^-^ 



=1 fcez" VmeZ" 



=1 fcez" \meZ' 

n n+l 



+ E E ( E '-^iU^{p')e-- 1 

= E 

feez" 

+ E f E f.'-^m^{p')e-A 

feez" Vmez" ;=i ^^^"+1 / 



EE( E ^m,^„,(iM|^(p')e^-^ 



3 = 1 1=1 \mg 



E 



71 71+1 



EE( E ^'^j^'"-iWf^J^^7^(p') 



7=1 1=1 VmeZ" 



+ E (e E^w^';|r'(^')l^""- 

wez" Vmez" i=i ^"+1 / 



n+l 

E 

j=l 



'"+1 QT- 



^i^(^')plE|f(p')e.We-- 



,i=i 



n+l 

E 



= E 

j=i 



E|.fE^MWe--)p;E|^(P')e.We-^ 

■"+1 ni, 
.1=1 fe6Z'» feez" 
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n+1 



=1 fceZ" 



= E 



ri+l 



1=1 feezi \mez" 



+ E E E ^(*Mj^(^^')e.We^-- e-- 



= E 



n n+1 / 

EE( E ^'^j^m,iit)p'i-Q^ip')ekit)e 



j—1 l—l 
^ n+1 



Jk-q 



+ E ( E E^W^^jJ^(^'')e.We--le*« 
fceZ" VmeZ" 1=1 ^Pn+i 



-{p')ew-m{t) 



n n+1 / Qj_^ 

E EE E ^"^j^^Mp'i^g^ 

E f E E^(*M^(^')e.-™w)e-- 

MiGZ" \meZ" 1=1 / 



4/? i/ie nineteen terms have one of three forms. Terms one, five, six, twelve, fourteen, and fifteen have the 
form 



E wy'"- 



Terms two, three , four, seven, eight, nine, ten, thirteen, sixteen, seventeen, eighteen and nineteen have the 
form 



J2 9k{p')hkit)e"' '. 



Term, eleven is simply a function of p' , say L{p'). Adding all terms together we see the final result, {x,F}, 
is of p' — q' — exponential-form. Finally we have have the following analyticity domains 



dr do dX 

Yj{q')€Ap',a', g^{p',q') € Ap-S,a, ^{p' ,q') & Ap,a-5, -g^{q')&Ap,,a'-5, 

1=1 ^3 Pj fcez" Pj 

for the functions that make up the expression for {XtF}{p\q')- Since these functions are multiplied together 
in a number of ways and then summed it follows, given p' < p and a' < a, {x, F}{p', q') S Ap'-s,a'-s o-nd is 
of p' — q' — exponential form. □ 
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Lemma F.8 

Let H[p' , q') = U{p' , q') + P{p' , q') with H and P beAng of p' , q' -exponential form. Define TZa = '''2{U, Xi 1) 
ri(P, X, 1), where 

m — 1 ,/ oo ,1 

rUH,x,t)=UH = E TT^x^' 

L'^H — H and L™i7 = {L™^^i/, H} for m> 1. TZa is of p' , q' exponential form. 
Proof: 



We have 



where 



T^A ^{x.P}+UH-H- {x, H} = r2(C/, x, 1) + x, 1) 



,iH,x,t)=UH^ E 7[^x^ = E 7t4^- 



7^A = {x,n + ^2(i^,x,l)• 

By lemma \F.1\ with x G ylp',o-' /or some positive 6 < a' it follows {x, P} G Ap'-s.a'-s is of p' — 
q' exponential form. Next we examine 

oo ^ 

r2(i?,x,l) = Ez!4^- 

1=2 

Clearly by inductively applying lemma [FV^ r2{H, X; 1) o/p' — q' exponential form in some domain, "Dp^ gj^, 
to be determine. We set 

oc 
1=0 



a* 

i=0 



Each time we apply lemma F.7 , we can choose Si arbitrarily small. In particular we can set Si = Si-xjl for 

oo /I \ » 



i — 1,2, ... with So ~ S. Therefore 



i=0 ^ ' 
oo , -.i 

i=0 ^ ^ 

1/Fe finally have TZa G -4p»^cr* is of p' — q' exponential form. □ 



G Rossby Wave Flow 

The Rossby wave flow is generated by a Hamiltonian of the form 

H{X,Y,t) = H°{X,Y,t)+eH^{X,Y,t) 

= A sin ko{X - cot) sin/oF + eH'^{X, Y, t). 
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where e > 0, A is the maximum velocity in the y-direction, (/cq, ^o) is the wave number vector, and cq is the 
phase speed of the primary wave in the x-direction. In a reference frame moving with the primary wave the 
transformation 

X = X - Cot, y = 

yields the Hamiltonian 

H{x,y,t) = H"{x,y,t)+eH\x,y,t) 

= —Coy + ^sin kox sin loy + eH^{x, y, t). 
In this frame of reference the vector field generated by the Hamiltonian H{x, y, t) has the following form: 

X = Co — Alo sin kox cos loy — s— — (x, y, t), 

dy 

y = Ako cos kox sin loy + e— — (x, y, t). 

ax 

We will consider the perturbation H^{x, y, t) — xy sech^at where a is the decay rate. 

The near integrable form of this problem makes it a candidate to apply the theorem presented in this 
paper provided the conditions of the hypothesis are satisfied. The first step in applying the theorem is to 
transform the Hamiltonian to action-angle variables. Clearly, by the Liouville- Arnold Theorem |^ one can 
construct locally for the one degree of freedom system generated by the Hamiltonian H°{x,y) a symplectic 
coordinate transformation to action-angle variables. The integrable system is characterized by two hetero- 
clinic connections between saddle-type equilibrium points on the curves y = 0, y = tt respectively and by the 
symmetry trajectory starting at (0, 7r/2) shown in Figure 1. These three structures divide the phase space 
in four regions. Action-angle variables can be found locally for each of these four regions as we will now 
describe. The action variable p is defined as 



P=^ f xdy, (G.l) 



which is the area of the region inside the level set with Hamiltonian H divided by 27r. Figure 2 shows the 
actions for the system in Figure 1 starting from the right at p = for the elliptic equilibrium point in region 
/ and rising to a maximum at H — — 7r/4 for the symmetry trajectory. Similarly for regions /// and IV 
the action falls from a maximum at the symmetry trajectory to zero for the equilibrium point in region IV. 
Note I{H) is invertible in each of the intervals corresponding to regions /, II, III and IV. To define the 
angle 9{x,y), let L denote a straight curve emanating from the elliptic equilibrium point in region / to the 
elliptic equilibrium point in region IV. We denote solutions of the integrable vector field starting on L by 
{x{t, s), y{t, s)) where a::(0, s) = Xo{s) and y(0, s) = yo{s) so for any point (x, y) on the orbit {x{t, s), y{t, s)), 
t = t{x,y) is the time it takes for the solution starting at {xo{s),yo{s)) to reach {x,y). Given T{II) is the 
period of the periodic orbit with constant H, we define the angle variable, 9{x,y) as 

e{x,y)^2n'-^, (G.2) 

where {x,y) Cz H ~ constant. Clearly by this definition action-angle variables can not be defined on the 
heteroclinic connections. 

Next we show the Hamiltonian in action-angle variables is analytic. 

Lemma G.l The transformation to action-angle variables is analytic. 

Proof G.l We begin by transforming the integral part of the Hamiltonian II^(x, y) in region I . In region I 
([Clp has the following form 



= - / x{y,H°)dy = - / arcsin 



y 



sin y 2 



dy, (G.3) 
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X 

Figure 1: Phase space of the vector field generated by the Hamiltonian H = —coy + Asmkoxsinloy with 
A = ko = Iq = 1.0 and cq — 0.5. 



p 




-1.5 -1 -0.5 

Figure 2: Action map for each of the four regions starting from the right /, //, ///, IV. 
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0.2 



0.5 



Figure 3: Frequency Map. 



where the tt/2 term translates the integral to the y-axis and j/max; J/min cire the values of y where the level 
set intersects the y axis. We consider a complex extension of the variables , He = H'^ + iHjm, and y, 
yc = y + iyim, to complex strips and for sake of simplicity do not specify their width. The integral in ( [G-Sl ) 
becomes a contour integral along the real axis with end points ymax> J/min- Let y — {YmimYmax) denote 
the interval of y values in region I. Since T{y,H^) and T^^iy, H'^) are continuous for He in a complex 
extension, D, of H^(y) and T{y,H'^) and T^oiTj, H^) are continuous for yc on the contour consisting of 
the real interval (j/max, ymin)? it follows p{H^) is analytic in D . Furthermore, by the inverse function 
theorem, except for the elliptic equilibria where p = 0. H^{p) is analytic on some strip of p. Similarly, the 
proof for regions II, III, IV follows. 



For the angle variable in region I (CS) becomes 



e{x,y) 



27r 



dy 



sin^ y - {H° 



y\2 

2> 



T{H) 



dy 



By the same arguments as for the action transformation and defining a complex extension ofx, xc = x+ixim, 
it follows 9{x,y) is analytic for some strip of xc and yc- 

When transforming the perturbation H^(x, y, t) — xy sech ^at to the action-angle variables found for H'^ , 
it is clear that the time dependent term sech ^at is the same after the transformation. By the Liouville- Arnold 
Theorem the action-angle transformation {x,y) —> (p{H^),9(x,y)) is invertible, (j),9) — )■ (x{p,6),y(ji,9)), and 
by the inverse function theorem x{p, 9) and y{p, 9) are analytic. □ 

Note, since the time dependence of the perturbation is sech ^at, the perturbation is of exponential order 
with respect to time as required by the hypothesis of the theorem. 

Finally the nondegeneracy condition given by det{d^ H'~' / dp^) — det{duj / dp) ^ must be satisfied. We 
numerically compute the period and plot the frequency as a function of the action giving the frequency map 
shown in Figure 3. This graph gives the frequencies for the closed orbits in regions /, // as well as the closed 
orbits in regions ///, IV. The frequency at p = corresponds to both elliptic equilibrium points in regions 
/ and IV. The frequency falls to zero for a p value corresponding to both heteroclinic connections as one 
would expect. The frequencies to the right of this value of p correspond to the closed orbits in regions // and 
III reaching a maximum for the symmetry trajectory. The nondegeneracy condition is therefore satisfied for 
all tori except for the symmetry trajectory. 

H Aside Calculations 



We take an aside to check (9.2) reduces to the appropriate expression if the time dependence of g/c(s) is not 
aperiodic but is indeed periodic with frequency ujq. Without loss of generality we assume ujq to be positive. 
In this case we write the Fourier series expansion oi gk{s) 



9k{s) 
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and substitute this in (p.2|) which results in 



Now let X]/gz si — Ihe any series whose sum is one. Then we have 



-i(\-k)t 



Si 



E 



^9M / e''"°''e*(^-'^)(''-*)ds 



^i{\-k-]-ujQl)t 



E 



i{\ ■ k + (jJoO 



1 - e 



E 



i{\ • k + woO 



-i{X-k~\'UJt)l)t 



9kl 



i{X ■ k + luqI) 



Note, that except for the time dependence inside the square bracket, [ H.l ) is almost in the form of the 
Fourier series for fk[t) 



(H.l) 



ilujQt 



(H.2) 



Comparing (H.l) and (H.2) we see that we can eliminate the time dependence in the square bracket by setting 



gki 



t(\ ■ k + loql) 

By summing on both sides of the equation, this is equivalent to 



E^(o)-' = E 
/fe(o)E^' = E 



Qkl 



lei lei 



i(\ ■ k + ljqI) 



In which case, (ff.l) reduces to 



/fc(o)-E 



lez 



i{X ■ k + ojqI) 



(H.3) 



(H.4) 
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with 



fkl 



i{\ ■ k + ljqI) 



Note that [HA] reduces to the expression for fk{Q) obtained in (-ff.3) 



We now must obtain an estimate for \fk{t)\. Since G{q,t), F{q,t) e Ap, the Fourier coefficients fk{t) and 
gk{t) must be exponentially small with respect to the index k. First we analyze the estimate for gk{t)- 



We have 



G{q,t) = E^feWe*^ 



We calculate \gki\ 



J2 T.dkie''^"*^'''"'- 



1 1 /■•^ r"^ r^/^ 



gu 



(2t 



'-TT J-S/ 



where the integral is along the real axis. Lift the integral as follows 



gki 



fir fSl2 
-SI 



1 1 r r f'^' . fci . fc„ .1 

"SJ_^ J_^J_s/2 \ki\ \kn\ \l\ 



dqds 



gki 



(27r)"5' 



(^1"^ ki kn I 

/ G{qi - i—p, ...,qn - ijY^p,s - i — cr) 

-itJ-S/2 



\ki\ 



\kn 



and 



l5fcd<e-l''^l''e-l'l-°niG(g,t)||p,.. 



Now, since 



gkit)^Y.9kie''^°' 

is 



to calculate the estimate we need to calculate the following 



1/2 
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where t € C and t indicates the complex conjugate. Then 

|giiwo*| _ g-wo'(Ini t) 
where Im t is the imaginary part of t. Assume |Im t\ < a — S. Then 

-Woi(Im t) < \u>ol{lm t)\ < \l\\u)o\{cr - S). 



We then obtain 



\9k{t)\ < Y.\9ki\\e''"'"\ 

< e-l'=l^||G(5,t)||,,<.^e-l'l-°^ 



< 2e-l'=l''||G(g,i)llp.^Ee""°'' 
2e-l'=l^||G(g,t)||,,, ' ^ 



X — g— '^0'5 

We proceed similarly to obtain an estimate for fk{t). We know 



We obtain from the previous estimates 

\l\ujoS 



|AWI<e-l'=ll|G(g,i)llp,.E^ 



e 



A • fc + lu)o\ 

The diophantine condition on the frequencies gives an estimate for the denominator of this expression. First 
recall that for fc € Z" wo define \k\ = Y^^=i 1^*1 ^^'^ ll^ll = ^^Pj \ki\. Furthermore, define the vectors h = {k,l) 
and A = (A,a;o). Consider first the case h = {k,l) ^ 0. The diophantine condition is given by 

|A-/i| >r||/i||-("+i). 

For any K, Y,5'>0 recall the inequality 
It follows that 



101 



\A.h\ 



< r 



< r-^\h\ 



< r" 



n+ 1 
eS' 



n+l 



.\h\S' 



We choose 



so that 



Wo + 1 



l|G(g,t)||,,.r- 



< ||G(g,t)||,,<.r- 



< \\G{q,t)\UT- 



1 / (n+l)(wo + l) 



eoJoS 




(n+ l)(a;o H 




ea;oi5 




(n + l)(u;oH 


rl) 


eu)o5 



n+l 



e-l'=l(p-4r¥T)^e- 



n+l 



e ,fc|(p-<5) ^g-|i|wo5 

-|/c|(p-5) 



n+l 



1 _ g— woi5 



For the case h = {k,l) = we use the assumption G{q,t) = 0, where G{q,t) represents the average of the 
function over the variables q and t. In this case , we use the following Lemma to show /o(i) or /oo are 
bounded. 

Lemma H.l 

Given G{q,t) is periodic with respect to q and t and has Fourier coefficients gki, goo = iff G{q,t) = 0. 
Proof: 



We have 



and 



G{q,t) 



2-K S 



' iez 



V2 



-s/2 



E Y9kie'^"''e''"^dtdq 



9kl 



ILOqI 



-1 s/2 




TT - 


- -s/2 




— TT- 



500 



□ 



Therefore since 



we have /oo = and given 



Tkl 



9ki 



i{X ■ k + lojo) 



fk{t) = J2 



9kie 



i{X ■ k + Iujq) 
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it follows 



i/oWl < E 



l5o;||e 



I € Z 
I ^ 



\luJo\ 



< \\G{q,t)\\p,^r 



1 f { n + l){uJo + 1) 



n+l 



1 - e-'^f'^ 



Now we want to check that, given G{q,t) is periodic with respect to time, (9.10) reduces to the appropriate 
expression for fk{t), (HA). Assuming G{q,t) is time periodic with frequency luq implies the Fourier coefB- 
cients /fc(i) are periodic functions with 



and 



9k{t) = ^9kit 



The Fourier transform of a periodic function is equals to a series of the following from 



(H.5) 



That is, the Fourier transform consists of equally spaced delta functions which are weighted according to the 
Fourier coefficients of the function. 



Substituting {IL5) in ( ^.10 ) gives the following 



fkit) 



00 + 2/3 



Gk (w) 



D+ip i{X - k + oj) 



oo+i/3 i{\-k + uj) 



, , „ - ( 



(H.6) 



We use the following property of the delta function for a function h{Lj) continuous ai uj — 



S{Lu)h{uj)duj = h{0) 



and ( |H.q ) becomes 



9kl giicJot 



i{X ■ k + Iujq) 
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We derive a bound for the sum ^^^mGZ" ^ l"*l^e "l^^ "^^h We examine the cases n = 1, n = 2 and 
derive an expression for the general case. 



n= 1 



We obtain two bounds, one for fc > and another for fc < 0. These bounds are equal and thus the sum 
Emez" e-l'"l''e-l'=-™l(''-2'5) is bounded for any k. 

Let k > and divide the range of m in three intervals, m<0, 0<m<k and m > k. We obtain 
bounds for the three intervals 

m<Q ^ ^ g-|m|pg-|fe-m|(p-25) ^ ^-\k\{p-2S) ^ g"!'"!'', 
m<0 m<0 



g-|fe|(p-25)g-|m|(25) 



0<m<fc { ^ g-|m|pg-|fe-H(p-25) ^ ^ 

0<m< 

< g-|fe|(p-25) ^ g-|m|(25)^ 



0<m</c 0<m<fc 



Therefore for > 



0<m 



g-|fe-m|(p-2<5)^ 



m>k k—m<0 



e-l™lp 



mG^ m<0 



^ g-|m|pg-|fe-m|(^,-25) ^ ^-\k\{p-25) 

m<C 

+ g-|fc|(p-25) ^ g-|m|(26) 
I 

k—m< 



0<m 

g-|fe-m|(p-25) 

fe-m<0 



TO>0 



Let A; < and divide the range of m in three intervals, m<k, k<m<0 and m > 0. We obtain bounds for 
the three intervals. 

m<k I ^ e"l'"l''e"l'="™l(''"^*^ < e"l'=l'' ^ ^-\k-m\ip-25) ^ 

m<k k—m>0 



g-|fe|(p-25) g-|m|(25) 



fe<m<0 fc<m<0 

E 

fe<m<0 

^ -|™l(2«) 

m<0 
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< e-l'^K''-^^) ^ e- 



TO > I ^~\m\p^-\k-m\{p-2S) ^ g-|fc|(p-25) ^ 



m>0 m>0 

Therefore for < 

-m|(p-2<5) 



mG^ fc— m>0 



e i-^l(2"5) 



^ g-|fc|(p-25) ^ 
m<0 

+ g-|fc|(p-25) ^ 
m>0 



-|fc|(p-25) \^ g-|m|p 
m>0 



m>0 



Finally for any k G Z we have the bound 



n = 2 



In this case we have vectors k — (fci,fc2) £ and m = {mi, 1112) E I?'. We must keep in mind four 
sub-cases (1) fci > and k2 > 0, (2) fci > and k^ < 0, (3) A;i < and fca > 0, (3) ki < and fcz < 0. We 
will find a bound for sub-case (1) and the rest will follow without loss of generality. We now find the bound 
for ki,k2 > 0. 



mi < 
m2 < 



TUj^ < TTlj^ < 

7712 ^ *-* "^2 ^ ^ 



\m\p 



™1 ^ / g-|m|Pg-|fe-m|(p-25) 



< m2 < 



E 

»i < 

< ^ g-|mi|pg-|fei|(p-25)g-|fe2|(p-25)g-|m2|(25) 



mi < 
< m2 < ^2 



mi < 
< m2 < ^2 



g-|fe|(p-25) g-|mi|pg-|m2|(25) 



mi < 
< m2 < k2 



g-|m|(25) 



< mi 
< m2 



mi<0 ( ^_|^|p^_|fc_^|(p_25) 
7712 > k2 I 

mi < 
m2 > fc2 
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< ^ g-|'ni|Pg-|fei|(p-25)g-|fe2|Pg-|fe2-m2|(p-25) 

mi < 
m2 > fc2 

< g-|fe|(p-25) ^ ^_|TOi|p^_|fe2-m2|(p-25)^ 

mi < 
m2 > k2 



< mi < ki 

7112 < 



{ E 



g-|"»|Pg-|fc-™l(p-2*) 



< e 



< e 



< mi < fci 
m2 < 

-\k\{p-26) 



E 



g-|m2|Pg-|mi|(25) 



< mi < fci 
m2 < 



-\k\{p-2S) 



E 

m2 < 



-\m2\p 



< mi < fci 
< 7712 < fc2 



< 



E 

: mi < 
: m2 < 

E 



< mi < fci 
< m2 < fc2 



g-\m\p^-\k-m\{p-2S) 



g-|fci|(p-25)g-|mi|(25)g-|fc2|(p-25)g-|m2|(25) 



< mi < fci 
< m2 < fc2 



= e 



-\k\{p-2S) 



E 



g-|mi|(25)g-|m2|(25) 



< mi < fci 
< m2 < fc2 



< e-|fc|(P-25) ^ g- 



|mi|(25)g-|m2|(25) 



< mi 
< m2 



< mi < fci 

m2 > fc2 



< 



E 

mi < 
,2 > f=: 

E 



< mi < fci 
m2 > fc2 



g-|"i|Pg-|fe-m|(p-2(5) 



g-|fci|(p-25)g-|mi|(25)g-|fc2|Pg-|fe-m2|(/9-25) 



< mi < fci 
m2 > fc2 



< e 



-\k\{p-2S) 



E 



g-|mi|(25)g-|fe2-m2|(/9-25) 



< mi < fci 
m.2 > fc2 



mi > ki 
m2 < 



{ E 



g-|™|Pg-|fe-™l(p-2'5) 



mi > fci 
m2 < 



< 



mi > fci 
m2 < 



< g-|fe|(p-25) 



g-|m2|pg-|fei-mi|(p-25) 



mi > fci 
m2 < 
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mi>ki J- ^_\^ip^_ik-m\{p-25) 



< m2 < I ^ 

mi > ki 
< m2 < ^2 



< ^ g-|fci|pg-|fc2|(p-25)g-|m2|(25)g-|fei-mi|(p-25) 



mi > fci 
< >7i2 < fc2 



< g-|fc|(p-2,$) ^ 



g-|m2l(25)g-|fci-mil(p-25) 



mi > fci 

< 7712 ^ ''-2 



"^1 ^ ^1 ) \ " g-|m|pg-|fe-m|(p-25) 



{ E 



1712 > ^2 

m2 > fc2 



< ^ g-|fei|pg-|fei-mi|(p-25)g-|fe2|Pg-|fe2-m2|(p-25) 



mi > fei 

m2 > fc2 



m|(p-25) 



mi > fci 
m.2 > fc2 



Finally we have 



mi < 

m2 < 

_^ g-|fc|(p-25) ^ g-|m|(25) 



mi G Z mi < 

m2 G Z m2 < 



< mi 
< m2 



\k\(p-2S) ST g-|mi|pg-|fe2-m2|(p-25) 



g-|'"2|p 



"2 < 



\k\(p-2S) Sr f,-\mi\(2S}^-\m2\{2S) 

< mi 
< ma 



^ g-|/c|(p-25) ^ 

mi < 

m2 > k2 

+ g-|fe|(p-25) 

m2 < 

^ g-|/c|(p-25) ^ 

< m 
< m 

+ g-|fe|(p-25) 

< m] 
m2 ^ 

+ g-|fe|(p-25) ^ 

mi > i 
m2 < 

+ g-|fc|(p-25) Y 

mi > ? 
< m2 < 



g-|mi|(25)g-|fc2-m2|(p-25) 



< mi < fei 
m2 > fc2 



|fe|(p-25) ^_lm2\p^-\ki-mi\{p-2S) 

mi > fci 
m2 < 

|fe|(p-25) g_|TO2|(25)g-|fei-mi|(p-25) 

mi > fci 
< m2 < fc2 



-2 > fc2 



|mi|(25)„-|m2|(25) 



I mi 
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Without loss of generality the same bound is obtained for the cases involving different values of ki and 
A;2 -Therefore we have for any k GiP 

g-|m|pg-lfc-m|(p-25) ^ g^-\k\{p-2S) f ^ \ 

Following the same procedure we can write a bound for the arbitrary n case. That is, for any A; e Z" 

^ g-|m|pg-|fe-m|(p-25) ^ ^n^-\k\{p-2S) f ^ 
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